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FREPACE 

Nonmathematlcal  Summary 


This  study  Is  concerned  with  the  problem  of  choosing  the  amounts 
x^jXg, . .  of  a  single  product  (or  aggregation  of  several  products) 

to  produce  In  each  of  n  successive  time  periods  l,S,,,,,n  so  as  to 
minimize  the  total  manufacturing  costs  over  the  n  periods.  The  re¬ 
quirements  ^i>^2''*‘'^n  product  occurring  in  periods  1,2, ,..,n 

are  assumed  to  he  known  in  advance.  Requirements  in  a  period  are 
satisfied  in  so  far  as  possible  from  stock  on  hand  at  the  beginning  of 
the  period  and  from  production  during  the  period,  Requirements  which 
cannot  be  met  in  this  way  are  backlogged  until  they  can  be  satisfied 
by  subsequent  production. 


Let 


J. 

yi  “  2  ^1  ^  h 

j  “1 


the  amount  of 


Inventory  on  hand  at  the  end  of  period  1.  If  y^^  is  negative,  then 
-y^^  is  the  total  amount  by  which  the  cumulative  requirements  exceed 
the  cumulative  production  in  the  first  1  periods.  We  suppose  that 
this  excess,  -y^^,  is  backlogged  until  it  can  subsequently  be  satisfied. 

Denote  the  cost  of  producing  x^  units  in  period  1  by  o^(x^), 
The  cost  of  storing  y^^C  >0)  units  at  the  end  of  period  1  is  denoted 
ty  h^(y^).  When  y^  <  0,  is  the  penalty  cost  Incurred  because 

-y^  units  of  requirements  are  backlogged  at  the  end  of  period  1. 

We  suppose  that  there  are  given  upper  and  lower  limits  x^^  and  x^ 
(x^  <  x^)  respectively  on  production  in  period  1  (=  l,2,...,n).  In 
addition  there  are  given  upper  and  lower  limits  y^^  and  y^  <  y^) 
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on  in  period  i  (=1,2, .. ,,n-l)j  and  no  inventory  or  ‘backlog  is 
allowed  at  the  end  of  period  n,  i.e.,  yj^  “  0. 

The  problem  is  to  choose  production  levels  x^, that 
minimize  the  total  coat 


S  n-1 

1  +  I 

i=i  ^  ^  1=1  ^  ^ 


over  the  n  periods  subject  to  the  above  named  constraints. 

We  assume  that  the  cost  fimctlons  c^  and  h^  are  convex.  For 
example  they  might  appear  as  in  Figure  1. 


Figure  1 


(The  functions  in  Figure  1  are  convex  because  the  chord  connecting  any 
two  points  on  the  graph  of  either  function  does  not  fall  below  the 
graph  of  the  function  between  the  two  given  points . ) 

The  objective  of  our  study  is  to  determine  the  effect  of  changes  in 
the  requirements  and  capacity  limitations  (production,  storage,  backlog) 
upon  the  optimal  production  levels.  We  show  that  the  optimal  production 
level  in  a  given  period  is  a  non-decreasing  function  of  (l)  the  require¬ 
ments  in  any  period,  (2)  the  upper  and  lower  production  capacity  limits 


ii 


in  the  given  period,  and  (3)  the  upper  and  lower  storage  limits  ^  in 
the  given  and  all  succeeding  periods;  the  optimal  production  level  in 
a  given  period  is  a  non- Increasing  function  of  (l)  the  upper  and  lower 
production  capacity  limits  in  any  other  period  and  (2)  the  upper  and 
lower  storage  limits  in  all  preceding  periods.  The  above  results  also 
lead  to  simple  and  efficient  computational  procedures  for  finding 
optimal  production  quantities. 

As  an  illustration  of  these  methods,  suppose  that  the  production 
levels  are  assumed  to  be  integers.  Suppose  also  that  we  have  optimal 
(Intogsr)  production  levels  ...,x^  for  the  (integer)  requirements 
^l'“*'^n*  optimal  production  levels  x^,  ...,x^  for  the  re¬ 
quirements  where  1  for  some  k  and  r|  =  r^  for 

1  /  k.  In  this  circumstance  we  show  that  there  is  an  integer  J  (not 
necessarily  equal  to  k)  for  which  x^  «  Xj  +  1  and  x^  =  x^  for  1  J. 

Example; 

As  an  example  of  the  way  in  which  this  result  can  be  applied  to 
find  an  optimal  production  schedule,  suppose  the  data  are  as  given  in 
table  1  below. 


ll 

In  this  discussion  backlogged  requirements  are  viewed  as  negative 
inventories. 
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Table  1 


The  idea  Is  to  begin  by  defining  a  sequoiiod  of  requirements  for  which 
there  is  only  one  feasible  set  of  production  levels.  For  example^  if 
■  -2,  Tg  °  0,  rj  =>  1,  rjj^  =  1  then  Xj  "  0  for  1b1,2,3,  4  is  the 
only  feasible  set  of  production  levels  since  «  0  for  all  i.  We 

now  increase  the  requirements  one  unit  at  a  time  until  we  obtain  the 
requirements  in  table  1.  At  each  stage  we  find  a  corresponding  optimal 
collection  of  production  levels. 

We  begin  by  increasing  the  requirements  in  (say)  period  3  to  two. 
From  what  we  have  said  above,  this  requires  us  to  produce  one  additional 
unit  in  one  of  the  first  four  periods.  Increasing  production  in  periods 
one  or  four  is  not  feasible,  the  former  because  the  upper  Inventory 
limit  in  period  one  would  be  violated  and  the  latter  because  the  upper 
limit  on  production  in  period  four  would  be  violated.  The  cost  of 
producing  the  one  unit  in  periods  two  and  three  is 
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eg(l)  +  h^(2)  +  hg(3)  +  hj(l)  =  13+3+8+1  =  25 

and 

cjd)  +  h3^(2)  +  hg(2)  +  hjCl)  =  9+3+5+1  =  18 

respectively.  Thus>  the  production  levels  =  1  and  x^  =  0,  1  ^  3> 

are  optimal  for  the  requirements  r^  =  -2,  =0,  r^  =  2,  r^^  =  1. 

Next  we  Increase  the  reqxiirements  In  period  two  hy  one.  This  time  the 
most  economical  plan  Is  to  Increase  production  In  period  two  by  one. 
The  final  step  Is  to  Increase  requirements  In  period  three  again  by 
one  unit  obtaining  -2,  rg  =  1,  r^  =  3>  r^^  =  1.  The  best  plan  is 
now  to  Increase  production  In  period  two  by  one,  obtaining  x^  =  0, 

Xg  =  2,  Xj  =  1,  Xj^  =  0  as  the  final  set  of  optimal  production  levels. 


HRODUCflON  PLANNING  WlOffi  CONVEX  COSTS: 

A  PARAMETRIC  STUDY 

by 

Arthur  P.  Veinott,  Jr. 

1.  Introduction  and  Summary 

We  Goneider  the  problem  of  choosing  the  amounts  of 

a  single  product  to  produce  In  each  of  n  successive  time  periods 
l)2>...,n  so  as  to  minimize  the  total  manufacturing  costs  over  the 
n  periods.  The  requirements  product  occurring 

in  periods  1,2, ...,n  are  known  in  advance.  Requirements  in  a  period 
are  satisfied  in  so  far  as  possible  from  stock  on  hand  at  the  beginning 
of  the  period  and  from  production  during  the  period.  Requirements  which 
cannot  be  met  in  this  way  (because,  for  example,  of  limited  production 
capacity)  are  backlogged  until  they  can  be  satisfied  by  subsequent 
production.  It  is  convenient  in  the  sequel  to  view  backlogged  require¬ 
ments  as  a  negative  inventory.  Similarly,  disposal  of  ("excess")  stock 
is  viewed  as  negative  production. 

We  admit  two  types  of  costs  in  a  given  period:  production  and 
holding,  the  former  being  a  disposal  cost  when  the  production  level  is 
negative  and  the  latter  being  a  backlogging  cost  when  the  inventory 
level  is  negative.  These  costs  are  assumed  to  be  convex  functions 
respectively  of  the  qiiantities  produced  during  and  stored  at  the  end  of 
the  period.  In  addition  we  permit  upper  and  lower  limits  to  be  imposed 
on  the  amovints  produced  and  stored.  The  cost  functions  and  quantity 
limitations  for  successive  periods  need  not  be  the  same. 


The  model  can  be  interpreted  in  a  variety  of  different  ways  including 
service  scheduling,  warehousing  decisions,  and  distribution  of  effort. 
Several  of  these  possibilities  are  discussed  in  [4]  and  [?]• 

Our  objective  is  to  study  the  effect  upon  the  optimal  production 
levels  of  changes  in  the  parameters,  i.e.,  the  requirements  and 
quantity  limitations.  In  the  event  that  the  total  cost  function  is 
strictly  convex,  the  optimal  production  quantities  are  unique.  For 
this  case  our  two  main  results  are  easily  stated  as  follows.  First, 
the  optimal  production  level  in  a  given  period  is  a  non-decreasing 
function  of  (l)  the  requirements  in  any  period,  (2)  the  upper  and 
lower  production  capacity  limits  in  the  given  period,  and  (j)  the 
upper  and  lower  storage  limits  in  the  given  and  all  succeeding 
perlodsj  the  optimal  production  level  in  a  given  period  is  a  non- 
increasing  function  of  (l)  the  upper  and  lower  production  limits  in 
every  other  period  and  (2)  the  upper  and  lower  storage  limits  in  any 
preceding  period.  Second,  the  optimal  cumulative  production  levels 
X^,X2, . . .,X^^  are  each  non-decreasing  functions  of  the  cumulative 


requirements  R^,R2>...^E^  in  each  period,  where  X^ 
1 


1 

y  X.  and 
J=1 


R.  =  ^  r.  for  all  1. 

^  j=l  ^ 

The  first  result  Is  of  interest  for  studies  of  possible  changes 
in  the  production  or  storage  capacities,  or  changes  in  the  minimal 
guaranteed  production  level.  Itie  first  result  is  also  the  basis  for 
some  efficient  parametric  programming  procedures  to  be  described  later. 
The  second  result  is  useful  in  forecasting.  For  example  suppose  that 
we  do  not  know  the  actual  cumulative  requirements  R^,R2, . . .,R^,  but 
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can  forecast  maximal  and  minimal  cumulative  requirements  R^,Rg, . . 
and  R^jRg, respectively  with  assurance  that  <,Rj  <  for 
all  i.  We  may  then  compute  the  corresponding  optimal  maximal  and 
minimal  cumulative  production  levels  X^,X2, and  X^^Xg, ...,2^ 
and  he  assured  that  the  cumulative  production  levels  X^,Xg, ...,X^  that 
are  optimal  for  Rj^,Rg, . .  .,R^  ,  are  such  that  for  all  1. 

At  the  beginning  of  period  one  we  will  ordinarily  only  have  to 
choose  and  not  subsequent  production  levels.  If  it  happens  that 
X^  =  Xj^  —  and  this  can  occur  even  when  the  and  R^  are  not  all 
the  same  —  then  the  optimal  amount  to  produce  in  period  one  is  determined 
as  X^  without  complete  knowledge  of  the  E^.  If  instead  X^  <  X^, 

then  Improved  forecasting  is  needed  in  order  to  determine  the  optimal 
value  of  X^.  The  value  of  narrowing  the  forecast  Interval  in  various 
periods  can  be  assessed  by  determining  its  effect  upon  the  width  of  the 
interval  [X^,  X^]. 

The  rather  intuitive  relations  described  above  between  the  optimal 
production  levels  and  the  various  parameters  do  not  hold  in  general, 
but  do  seem  to  be  valid  in  a  number  of ■ situations  not  encompassed  by 
our  hypotheses.  It  is  somewhat  surprising,  however,  that  the  results 
may  fail  to  hold  when  the  total  costs  are  convex  but  cannot  be  expres¬ 
sed  in  the  form  assumed  in  this  paper.  For  example,  suppose  the  total 
cost  for  a  two  period  model  is 

max(xg,5x^  +  2x2-6)  . 

This  fuinctlon  is  convex  since  the  maximum  of  two  linear  functions  is 
convex.  However,  if  the  requirements  are  r^=  r2=  1,  then  the  optimal 


5 


production  levels  are  x^=  2,  Xg  =  0,  while  if  r^=  1,  rg=  3  the 
optimal  production  levels  are  Xj^=  1,  Xg=  3  (assuming  Xg  =  Rg  and 
x^  >  0,  Xg  >  0).  Observe  that  neither  of  our  two  results  holds  here. 

The  fact  that  the  optimal  production  quantities  are  non-decreasing 
functions  of  the  requirements  leads  Immediately  to  an  extremely  simple 
procedure  for  computing  optimal  production  schedules.  In  order  to 
Justify  the  algorithm  we  are  about  to  describe,  it  is  necessary  to 
make  two  additional  assumptions.  First,  the  cost  functions  are  piece- 
wise  linear  with  the  endpoints  of  each  of  the  intervals  on  which  the 
functions  have  linear  segments  being  Integers  also.  Second, the  require¬ 
ments  and  the  upper  and  lower  limits  on  production  and  storage  are  all 
integers,  Under  these  conditions  we  show  that  there  are  optimal 
production  levels  that  are  also  integers.  ^ 


In  the  most  elemental  problem  with  which  our  algorithm  deals,  we 
start  with  optimal  production  levels  x^^,  ...,x^  (all  integers)  for  a 
given  sequence  of  requirements  r^, We  then  wish  to  find  optimal 
production  quantities  x£, ...,x^  for  the  new  requirements  r^, ...,r^ 
where  r^  =  r^  +  1  (=  r^-l)  for  some  1  and  r^  =  r^  for  J  ^  1. 

Our  basic  result  is  that  there  is  an  integer  k(l  <  k  <  n)  for  which 


x^  =  Xj^  +  1  (=  Xj^-1)  and  X- 


X  for  J  k. 


This  means  that  in 


This  fact  is  not  surprising  since  several  authors  [l],  [6]  have 
shown  that  special  cases  of  our  model  can  be  formulated  as 
transportation  type  linear  programming  problems  for  which  Integer 
solutions  exist.  (The  amounts  "shipped"  and  "received"  must  be 
Integers,  of  course,  as  they  would  be  under  our  assumptions.) 

Our  problem  can  also  be  formulated  in  this  way. 


k 


order  to  find  {x^  ,  It  suffices  to  compute  the  total  cost  under  each  of 
the  n  possible  production  sequences  (corresponding  to  k=l,2, n), 
and  then  to  choose  the  cheapest  one. 

I 

Our  approach  to  more  complex  problems  involves  solving  a  sequence 
of  elemental  problems.  For  example,  suppose  that  we  seek  optimal  pro¬ 
duction  levels  for  a  requirements  sequence  rj^,  ...,r^  for  , 

which  the  r^  are  arbitrary  integers.  We  begin  by  defining  a  sequence 
of  requirements  rj^,...,r^  for  which  there  is  only  one  feasible  set  of 
production  levels  x^, ...,x^.  To  illustrate,  if  disposal  of  stock  is  not 
allowed,  we  may  choose  x^=t^=0  for  all  J  (provided  that  this  choice  is  ^ 

feasible).  We  then  construct  a  new  requirements  sequence  by  adding 

+l(-l)  to  any  r,  for  which  (r'-r.)  is  positive  (negative).  Next  we  find 
J  |1  J 

an  optimal  production  schedule  for  the  new  requirements  sequence  using 

the  technique  described  in  the  preceding  paragraph.  We  then  repeat  the 

process  after  replacing  the  original  Tj  by  the  new  requirements.  We 

continue  in  this  way  until  we  find  that  rj  =  r^  for  all  J,  at  which 

point  the  original  problem  is  solved.  The  process  must  terminate  because 

n 

at  each  stage  we  reduce  |r'-r  |  by  +1. 

J=1  ^  ^ 

When  the  r!  are  all  non-negative,  one  can  view  the  above  process  as 
0 

satisfying  each  unit  of  requirements  in  turn  as  cheaply  as  possible.  The 

t 

special  case  of  this  procedure  in  which  requirements  must  be  satisfied  in 
order  of  their  due  dates  is  shown  by  Johnson  [2]  to  be  optimal  for  the 
special  case  of  our  problem  in  which  no  backlogging  is  allowed,  no  storage 
limits  are  permitted,  and  the  Inventory  carrying  costs  are  linear. 

The  parametric  pix)gramming  procedure  that  we  have  described  above  in¬ 
volves  only  changes  in  the  requirements.  In  section  4  we  develop  a  generalized 

the 

procedure  in  which  changes  in/productlon  and  storage  limits  are  allowed  as  well. 
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Alternative  procedures  for  solving  the  problem  considered  by- 
Johnson  are  offered  in  [It],  [5],  and  [8].  In  [9]  a  generalization  of 
Johnson's  procedure  is  shown  to  be  applicable  where  the  price  received 
for  each  unit  of  product  sold  depends  on  the  qioantitjr  sold.  ^  A  dynamic 
programming  procedure  that  solves  the  general  problem  in  this  paper  and 
that  neatly  exploits  the  convexity  assumptions  is  given  in  [4], 

References  to  earlier  work  will  be  found  in  [1],  [^t],  [5],  [6],  [7], 

[8],  [9]. 

We  outline  the  plan  of  this  paper.  The  problem  is  formulated  in 
section  2.  We  establish  certain  fundamental  inequalities  in  section  3. 

In  section  4  we  develop  computational  procedures  for  finding  optimal 
production  levels. 

2.  Formulation  of  the  Problem 

Let  y^  be  a  given  constant  and  let 
1 

(1)  ^i  '  S  ^o  ^  “  l,2,...,n  . 

We  interpret  | y^ |  as  the  total  amount  of  stock  on  hand  at  the  end 
of  period  1  when  y^  >  0,  and  as  the  total  amount  of  backlogged 
requirements  at  that  time  when  y^  <  0.  We  are  implicitly  assuming 


Actually  our  model  provides  for  this  possibility  by  allowing 
disposal  of  stock.  We  omit  a  discussion  of  this  point,  however. 
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that  it  is  not  possible  to  end  a  period  with  both  a  positive  inventory 
and  a  positive  amount  of  backlogged  requirements.  Similar  remarks  apply 
to  production  and  disposal  in  a  given  period. 

There  are  given  upper  and  lower  bounds,  y^  and  y^  (<  y^)  on 
i.e., 

(2)  5  5  i  =  1,2, , . ,,n-i  . 

Observe  that  by  choosing  y^^  =  0  for  all  i  that  we  eliminate  all 
backlogging.  We  may  also  allow  minimum  inventory  levels  (e.g.,  as  a 
hedge  against  uncertainty)  by  letting  y^  be  positive. 

We  suppose  that 

(3)  for  1=1,2, ...,n 

where  the  x^^  and  x^  (x^  <  x^)  are  given  constants.  The  presence  of 
minimal  production  levels  enables  us  to  study  the  effect  of  different 
guaranteed  employment  policies. 

Let  g(x,r)  be  the  total  cost  associated  with  the  production 
schedule  x  =  (x^)  and  requirements  schedule  r  =  (r^^).  ^  We  assxme 
that  g  can  be  written  in  the  form 

(^)  g(x,r)  =  2  i  (y^) 

i=l  ^  ^  i=l 


2/ 


'bll 

We  write  z  =  (zj^)  to  mean  that  z  is  a  vector  whose  i 
coordinate  is  . 
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where  and  are  convex  and  continuous  on  the  intervals 
[Xj^,  Xj]  and  [y^,  y^]  respectively.  In  applications  c^(x^)  is 
the  cost  in  period!  of  producing  x^  units  if  >  0  and  of  disposing  of 
-Xj^  units  if  x’^  <  0}  h^(y^)  is  the  cost  in  period  1  of  carrying  y^ 
units  in  inventory  if  y^^  >  0  and  of  backlogging  -y^  units  if 
<  0. 

There  la  no  loss  of  generality  in  assuming  that  y^  =  0  since> 
by  virtue  of  (l),  (2),  and  (4),  we  may  replace  r^  by  (r^  -  y^). 

This  substitution  could  lead  to  a  negative  value  of  r^^  -  a  possibility 
that  we  do  not  rule  out.  Indeed  we  permit  any  of  the  r^  to  be  negative. 

Our  final  assumption  is  that 


which  states  that  we  must  end  period  n  with  no  inventory  or  backlog. 

It  is  convenient  to  let  x  =  (x^)  ,  x  =  (x^)  ,  y  =  (y^)  , 

Z  -  P  P  ®  parameter  set. 

In  the  seguel  we  shall  say  that  the  production  schedule  x  is 
feasible  for  the  parameter  set  p  (or  briefly,  x  is  feasible  for 
p)  if  the  restrictions  (l),  (2),  (5),  and  (5)  are  satisfied.  If  x 
minimizes  (4)  subject  to  the  above  constraints,  we  say  that  x  is 
optimal  for  p.  Finally,  we  say  that  p  Is  feasible  if  there  is  an 
X  that  is  feasible  for  p. 

Since  the  collection  of  feasible  production  schedules  Is  compact 
whenever  p  Is  feasible,  and  since  g  is  continuous,  there  exists  an 
optimal  X  whenever  p  Is  feasible. 


8 


In  certain  situations  it  Is  desirable  to  drop  the  assumption  (5). 
V?hen  this  is  the  case  ve  shall  suppose  that  there  are  given  limits 
and  on  the  inventory  level  at  the  end  of  period  n.  We  can 
then  imbed  the  original  n-period  problem  In  an  (n+l) -period  problem 
in  which  (5)  does  hold  without  loss  of  optimality. 

In  the  extended  problem  we  put  yJ^^  *n+l“  ^n"  ^n’  ^n+1  “ 

and  let  ^  ^  ^ 

feasible  production  schedule  for  the  n-period  model  in  which  (5)  is 
not  Imposed.  Then  in  order  for  the  schedule  x  =  (xj  'to 

feasible  for  the  (n+l)-perlod  model  in  which  (5)  does  hold,  we  must  have 


'n+l 


n+l 


n+l 


by  (1).  Hence,  -  y^^.  This  uniquely  defined  value  of 

Xn+i  also  satisfies  (5)  for  i  =  n+l  by  virtue  of  the  definition  of 
the  parameters  for  period  n+l. 

We  now  show  that  x*  is  an  optimal  schedule  for  the  n-perlod 
model  in  which  (5)  heed  not  hold  if  and  only  if  x*  =  (x*,  is 

optimal  for  the  (n+l)-period  model  with  (5)  holding.  First  suppose 
that  X*  is  optimal  for  the  n-period  model  and  that  x  =  (x, 
is  any  feasible  schedule  for  the  (n+l) -period  model  with  requirements 
schedule  r  =  (r,  Then  x  is  feasible  for  the  n-period  model. 

Also  g(x,r)  =  g(x,r)  >  g(x*,r)  =  g(x*,r),  which  was  to  be  shown. 

Wow  suppose  that  x*  is  optimal  for  the  (n+l)-period  model.  Then 
g(x,r)  =  g(x,r)  >  g(x*,r)  =  g(x*,r),  which  completes  the  proof. 

The  above  remarks  enable  us  to  confine  our  discussion  in  the  sequel 
to  situations  in  which  (5)  holds. 
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5.  The  Fundamental  IneqvialltleB 


In  this  section  we  establish  two  Important  inequalities  that  relate 
the  difference  between  two  parameter  sets  to  the  difference  between  the 
corresponding  optimal  production  schedules.  In  order  to  state  the  first 
result  it  is  convenient  to  introduce  a  definition.  Let  z  =  (z^)  be 
an  (n-l)  coordinate  vector  of  real  numbers.  Let  m(z)  be  the  corre¬ 
sponding  n  coordinate  vector  of  real  numbers  whose  i'^^  (l  <  i  <  n) 
coordinate  is  max  z.  and  whose  n^h  coordinate  is  zero. 

JSl  J 

Theorem  1 

Suppose  that  p  =  (-r,  x,  x,  y,  and  p'  s  (-r*,  x',  xS  y' ,  ^' ) 
are  feasible  parameter  sets  and  that  P'  <  P»  ^  If  x(x')  is  optimal 
for  p(p')>  then  there  is  an  x'(x)  that  is  optimal  for  p'(p)  and 
that  satisfies 

(6)  (x'-x)  >  -(x-x' )-(x-x' )-in(y-y' )  . 

Before  proving  the  theorem  we  note  briefly  some  of  its  implications. 
Observe  that  if  p  and  p'  coincide  except  that  r'  >  r,  then  (6) 


2/ 


Let  V  =  (vj^)  and  w  =  (w^^)  be  vectors  of  real  numbers.  We  say 

that  V  >  w  if  V.  >  w.  for  all  1. 

—  1  —  i 
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states  that  x'  >  x.  On  the  other  hand  If  p  and  p'  are  identical 
except  that  xj^  <  x^^  for  some  i,  then  (6)  states  that  xj  >  Xj  for 
J  /  i  and  x|  >  x^-(x^-xp.  Finally,  if  p  and  p'  differ  only  in 
that  yj^  <  y^,  for  some  i,  then  (6)  states  that  x^  > 

J  <  i  and  xj  >  Xj  for  J  >  i. 

We  remark  that  the  ineijuality  (6)  can  he  sharpened  at  the  expense 
of  complicating  its  statement.  The  simpler  treatment  is  adequate  for 
our  purposes,  however. 


Proof; 

It  is  convenient  to  let  y  =  (y^)  and  z  =  x-(x-x' )“(x-x' )- 
m(y-y' )-m(y-y' )• 

We  prove  only  the  first  part  of  the  theorem,  i.e.,  if  x  is 
optimal  for  p,  etc.  (The  proof  of  the  second  part  is  similar  and 
will  not  he  detailed.)  We  do  so  hy  showing  that  for  any  schedule 
X*  that  is  feasible  for  p'  and  that  does  not  satisfy  x*  >  z,  there 
is  an  alternative  schedule  x'  such  that 


This  special  case  of  theorem  1  was  motivated  in  part  hy  theorem  2' , 
p.  2t5  in  [5],  which  applies  to  a  generalization  of  our  problem  in 
which  the  requirements  in  each  period  are  random  variables. 

Theorem  2'  in  [5]  is  not  as  strong  as  our  theorem  1  for  the  case 
of  deterministic  requirements.  In  our  notation  theorem  2'  asserts 
only  that  >  X]_.  In  theorem  2  of  this  paper  we  shall  show  that 
the  inequality  xj^  >  x  holds  under  the  weaker  condition  that 
R|  >  for  all  1. 
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(l)  x'  is  feasible  for  p'j 

(li)  g(x',r' )  <  g(x*,r' ); 

(lii)  (x'-x,x'-z,y'-y)  has  one  (or  more)  fewer  non-zero  coordinates 
than  (x*-x,x*-z,y*-y).  ^ 


Assuming  for  the  moment  the  truth  of  the  above  result,  observe 
that  either  x'  >  z  or,  upon  replacing  x*  by  x',  that  it  is 
possible  to  construct  a  new  x'  having  the  properties  (l),  (ii), 

(iii).  Since  at  each  repetition  of  this  procedure  we  either  terminate 
by  finding  an  x'  satisfying  (i)  and  (ii)  and  having  the  property  that 
x'  >  z,  or  we  reduce  the  number  of  non-zero  coordinates  of  (x'-x, 
x'-z,y'-y)  by  one,  we  must  terminate  in  at  most  3n-l  steps.  To  complete 
the  proof  of  the  first  assertion  of  the  theorem,  it  la  sufficient  to 
observe  that  if  x*  is  optimal  for  p',  the  x'  that  is  obtained  at 
the  termination  of  the  above  process  satisfies  the  conclusions  of  the 
theorem. 

It  remains  to  describe  a  procedure  for  constructing  a  schedule  x' 
with  the  properties  (l),  (ll),  (ill).  Since  x*  ^  z,  there  exists  an 
i,  1  <  i  <  n,  for  which  x*  <  z^.  It  is  convenient  now  to  consider 
two  cases,  y|_^  <  y^_j^  and  y*_^  >  y^_^  . 


ll 

In  order  to  conserve  space  we  have  not  stated  explicitly  that 
y*  is  the  vector  of  inventories  associated  with  x*  and  r' 
because  the  notation  makes  this  fact  evident.  For  the  same 
reason  we  shall  subsequently  refer  to  yt  without  identifying 
it  as  the  l^li  coordinate  of  y*.  We  shall  follow  a  similar 
practice  hereafter  when  no  ambiguity  will  result. 
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Case  l! 


Let 


xj>x 


j' 


j  be  the  smallest  Integer  greater  than  i 
Such  an  integer  must  exist  for  if  not 


for  which 


=  ^.1  + 


<  ^i-i  + 


t=i 


which  contradicts  (5).  The  Inequality  follows  from  the  fact  that 

yj.]_  5  5  -  ^t  ^  ~  i+2,  ...,n  > 

and  r^  >  r^  for  all  i. 

Let  e  =  mln(z^-x|,  x|-Xj).  It  follows  from  the  preceding  die 
cussion  that  €  >  0.  In  addition,  as  we  now  show, 


(7)  Xj^-x*  >  €}  y^-y|  >  €  fox  t»l,i+i, . , J-l}  and  xJ-Xj  >  e  . 


The  first  and  last  inequalities  follow  from  the  definition  of  e.  The 
intermediate  inequalities  follow  from  the  first  inequality  and 


(8)  y^-y*  =  -  S 


I 

k=i 


Now  define  an  alternative  production  schedule  x'  by 
x'  =  X*  +  e(u^-Uj),  where  u^^  is  the  k^*^  coordinate  unit  vector, 
i.e.,  the  k^^  coordinate  of  u^^  is  +1  and  all  other  coordinates 
are  zero.  We  show  that  x'  is  feasible  for  p'  as  follows.  Using 
(6),  (8),  the  definition  of  e,  the  hypotheses  of  the  theorem,  and 
the  feasibility  of  x  and  x*  respectively  for  p  and  p',  we  have 
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(9) 


Xj  <  xj  <  +  € 


and 


<  2fj  <  <  xj-c  =.  X-  <  xj  <  5-  . 

Thus  since  x|  ■  x*  for  t  /  i,J,  x*  satisfies  (3).  Also  for 
i  <  t  <  J  , 

^  yt  <  ^  ^  "•  "'l  "  ""i  ^  ^ 

<  +  x|  -  z^-(y^-yp  +  c 

<  X*  -  Zi  +  y^  +  e  <  y^  * 

For  all  other  t,  y^  =  y||,  so  that  x'  satisfies  (2)  and  (5). 

Thus  x'  satisfies  (i). 

In  the  following  we  denote  by  D'*'f(x)  and  D"f  (x)  respectively 
the  ri^t  and  left  hand  derivatives  of  a  function  f  at  the  point  x. 
When  the  function  is  convex,  as  will  be  the  case  in  the  succeeding 
discussion,  we  may  be  assured  that  the  right  and  left  hand  derivatives 
exist. 

Using  the  convexity  of  c^  and  h^  and  (7)  we  have 
g(x*,r')  -  g(x',r')  =  c^(x*)  -  c^(x|+e) 

j-i  .  . 

>  e[-D"Cj^(x*+e)  -  ^  ° 

k— i 

>  e[.D-c^(x^)  -  I>\(yj^)  +  D'"cj(xj)]  . 
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In  view  of  thlQ  Inequality  we  may  establish  property  (ll)  by 


showing  that 


To  prove  (10)  we  show  first  that  the  schedule  x(8)  =  x-8(u.  -  u.) 
is  feasible  for  p  provided  that  0  <  B  <  e#  Assuming  for  the 
moment  that  we  have  done  this,  we  note  that  x  Is  optimal  for  p 
and  therefore  that 

S(x(8),r)-g(x,r}  ^  q  0  <  B  <  e  . 


Letting  B  -* 0+  ,  the  above  inequality  becomes  (lO), 

It  remains  to  establish  the  feasibility  of  x(B).  Using  (6),  (8), 
(9)>  the  definition  of  e,  the  hypotheses  of  the  theorem,  and  the 
feasibility  of  x  and  x*  respectively  for  p  and  p' ,  we  have 


and 


<  x^-e  <  x^(B)  <  <  x^ 


ij  <  -J  <  =,(«)  <  =■,  +  e  <  x*  <  < ; 


Also  for  1  <  t  <  J  , 


'i  ^  *i  ^  5:t  ■  *^1  ^  ^'t "  *i 
z^  +  x|  <  yt  -  €  <  y^(8)  <  y^  <  y^  . 
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For  t  /  l,i,  x^(B)  =  and  for  t  <  i  and  t  >  J,  y^(6)  =  so 

that  x(8)  is  feasible  for  p  as  claimed. 

To  establish  (iii)  ve  observe  from  the  definition  of  e  that  all 

but  the  i^^  and  coordinates  of  (x'-x)  and  (x*-x),  and  of  (x'-z) 

and  (x*-z)  are  identical.  Also  all  but  the  1^^  through  the  (j-l)'*^'^ 

coordinates  of  (y'-y)  and  (y*-y)  are  the  same.  On  the  other  hand 

the  i*^  and  coordinates  of  (x*-x)  and  (x*-z),  and  the 
til 

through  the  (j-l)  coordinates  of  (y*-y)  are  all  non-zero.  But  at 

least  one  of  these  coordinates  in  (x'-x),  (x'-z),  or  (y'-y)  is  zero, 
which  proves  ( Hi ) . 


Case  2i  yj_^  >  y^_^ 

Let  J  be  the  largest  integer  less  than  i  for  which  x*  >  xy 
Such  an  integer  must  exist  for  otherwise 


1-1 


1-1 

I 

t=l 


(xj 


't  t 


i-1 

)<  I 

t=l 


^’'t'^t^  ^i-1 


if  i  >  2  and  ~  Vq  ~  ^  i  =  !•>  which  is  a  contradiction. 

Let  G  s  min(z^-x*,  ^cJ-Xj,  yj_2”^i  l^  ’  *^l®arly  g  >  0.  In 
addition 

(11)  ^  '  ^“^t  -  ^  J  <  t  <  1-1  >  ^  ■ 

The  first  and  last  IneqLuallties  follow  from  the  definition  of  e.  It 
remains  only  to  observe  that  for  j  <  t  <  i-1  , 


l6 


i-1 


i-1 


k=t+l 


which  establishes  (ll). 

We  define  an  alternative  schedule  x'  for  p'  by  x'  =  x*+e(u^-Uj)t 
The  proof  that  x'  satisfies  (3)  is  precisely  the  same  as  for  case  1. 

Now  for  J  <  t  <  1,  we  have  using  (ll),  the  hypotheses  of  the  theorem, 
and  the  feasibility  of  x  and  x*  respectively  for  p  and  p'  that 


-  ^t  -  ^t  -  “  ^t  ^t  -  ^t  > 

while  for  all  other  t,  y*  =  y^,  so  that  x'  satisfies  (2)  and  (5). 
Thus  property  (i)  Is  established. 

If  we  now  employ  (ll)  and  the  convexity  of  the  c^  and  h^,  we 
find  that 


g(x*,r‘)  -  g(x',r')  =  °j(x^)  ” 
i-1 

+  i-1  + 

>  e[D  Cj(x*-e)  +  X  ®  "  ^'°i(x*+e)] 

k=J 

.  i-1  , 

>  e[D  Cj(xj)  +  X  1*  • 

k=J 

In  order  to  verify  property  (ii)  it  therefore  suffices  to  show  that 


(12) 


1-1 


D^Cj(xj)  +  D\(y^)  -  D-c^(x^)  >  0 


Let  x(8)  ■=  x-8(u^-Uj).  We  show  that  x(8)  is  feasible  for  p 
provided  that  0  <  8  <  e.  The  fact  that  x(8)  satisfies  (5)  follows 
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from  the  same  argument  used  In  case  1.  For  J  <  t  <  1,  we  have  from 
(11),  the  hypotheses  of  the  theorem,  and  the  feasibility  of  x  and 
X*  respectively  for  p  and  p'  that 

-  ^t  -  ^t^®^  5  5  ' 

and  y^(6)  =  y^  for  all  other  t.  Ihus,  x(B)  is  feasible  for  p. 
Since  X  is  optimal  for  p, 

>0  for  0  <  8  <  e  . 

Letting  8  -» 0+  ,  we  obtain  (12). 

The  proof  that  x'  satisfies  (ill)  is  exactly  the  same  as  for 
ease  1  upon  interchanging  the  roles  of  1  and  j .  This  completes 
the  proof. 

In  order  to  state  and  prove  our  next  theorem  it  will  be  convenient 
to  Introduce  notation  for  the  emulative  production  and  requirements 
schedules  associated  respectively  with  the  production  and  requirements 
schedules  x  and  r.  In  particular  let 

X  (X^,X2, ...,  and  R  =  (R^,R2,  ...,  ^j^)  • 

Theorem  2 

Suppose  that  p  s  (-r,  x,  x,  y,  ^)  and  p'  =  (-r',  x,  x,  y,  ^) 
are  feasible  parameter  sets  and  that  R'  >  R.  If  x(x' )  is  optimal 
for  p(p' )/  then  there  is  an  x'(x)  that  is  optimal  for  p' (p)  and 
that  satisfies  X'  >  X. 
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Proof: 


The  proof  of  the  first  assertion  In  the  theorem  (l.e,,  If  x  Is 
optimal  for  p,  etc.)  consists  In  showing  that  for  any  schedule  x* 
that  la  feasible  for  p'  and  for  which  X*  ^  X,  there  Is  an  alternative 
schedule  x'  that  has  the  properties  (l),  (il)/  (ill)  (let  z  =  x 
In  (ill))  given  In  the  proof  of  theorem  1.  Once  this  construction  Is 
justified,  the  first  assertion  of  the  theorem  follows  from  an  obvious 
adaptation  of  the  corresponding  part  of  the  proof  of  theorem  1.  It 
remains,  therefore,  to  develop  a  procedure  for  constructing  a  schedule 
x'  with  the  properties  (l),  (ii),  (ill). 

Since  X*  ^  X,  there  is  an  Integer  1  for  which  XJ  <  X^.  Let 
1  be  the  smallest  such  integer.  Denote  by  J  the  smallest  Integer 
greater  than  1  for  which  X¥  >  X..  The  integer  J  exists  because 
by  (5) 


>  0  . 

n  n  n  n  — 


Let  e  =  min[x.-x*,  xf-x.,  min 
^  ^  'J  "J  i<k<  J 


(yj^-y*)].  We  have  s  >  0  since 


x^-x*  =  -  (XJ-X^^_^)  >  0 

^k'^k  =  ^  °  for  1  <  k  <  j  . 

Now  define  the  alternative  schedule  x'  for  p'  by  x'  =  x*+e(u^-Uj). 
Thus  y^  =  yj+e  for  k  =  1, i  H, . ..,  j-1  and  y^  =  y*  otherwise.  Using 
these  calculations  it  follows  easily  that  x'  satisfies  (i). 
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It  is  now  possible  to  show  that  (il)  and  (ill)  hold  using  precisely 
the  same  argument  as  that  used  to  establish  the  same  properties  in  case 
1  of  the  proof  of  theorem  1.  This  completes  the  proof  of  the  first 
assertion  of  the  theorem.  The  proof  of  the  second  assertion  of  the 
theorem  follows  in  a  similar  manner. 

!(•.  A  Barametric  Programming  Procedure 

Theorem  1  plays  a  central  role  in  this  section  in  the  development 
of  procedures  for  finding  optimal  production  schedules.  Our  techniques 
are  especially  useful  when  we  seek  optimal  production  schedules  for 
several  different  parameter  sets.  This  is  because  we  use  the  informa¬ 
tion  gained  in  determining  an  optimal  production  schedule  for  one 
parameter  set  in  an  efficient  way  to  reduce  the  computations  needed  to 
flnt'  an  optimal  production  schedule  for  a  different  parameter  set. 

In  order  to  simplify  the  exposition  in  this  section,  it  is  con¬ 
venient  to  impose  the  following  assumption: 

P:  and  h^(y^)  are  each  piecewise-llnear  functions 

with  the  endpoints  of  each  of  the  intervals  on  which  the 
functions  have  linear  segments  being  Integers  also. 

The  basic  result  on  which  the  computational  procedures  of  this 
section  rest  is 

Theorem  5 

Suppose  that  p  s  (-r,  x,  x,  y,  ;^)  and  p'  =  (-r',  x' ,  x' ,  y' ,  ) 

are  feasible  parameter  sets  with  Integral  coordinates,  that  (p-p' )  is 
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a  unit  vector,  that  P  holds,  and  that  x(x' )  has  Integral  coordinates 
and  is  optimal  for  p(p')* 

(a)  If  (y',2'')  =  (y>  y)>  then  there  is  an  x' (x)  that  is  optimal 
for  p'(p)  and  that  has  the  property  that  (x'-x)  +  (x-x')  +  (x-x') 
is  a  unit  vector. 

(h)  If  (y-y')  +  (y-y')  =  for  some  k,  then  there  is  an 

x'(x)  that  is  optimal  for  p' (p)  and  for  which  either  (x'-x)  = 

.  j  1 

for  some  i  <  k  <  or  (x'-x)  =  0. 

We  defer  the  proof  of  the  theorem  briefly  in  order  to  explore  some 
of  its  implications.  The  value  of  part  (a)  is  that  when  we  are  given 
an  X  that  is  optimal  for  p,  we  can  find  an  x'  that  is  optimal  for 
p'  (where  the  upper  and  lower  bounds  on  Inventory  levels  are  unchanged) 
by  considering  only  the  n  production  schedules  obtained  by  separately 
adding  each  of  the  n  unit  vectors  to  x-(x-x')  -  (x-x').  On  the  other 
hand  if  we  are  told  that  x'  is  optimal  for  p',  part  (a)  assures  us 
that  we  need  compare  only  the  n  production  schedules  obtained  by 
separately  subtracting  each  of  the  n  unit  vectors  from  x'  +  (x-x')  + 
(x-x' ) 


Example  1: 


If  X  =  (5,0,l)  is  optimal  for  p  where  r  =  (-2,6,0),  x  =  (3,3,5), 

X  =  0,  y  =  (5,7),  and  ^  =  (0,-2),  and  if  p'  is  such  that  r'  =  (-1,6,0) 
and  (x',  x',  V' >  ^')  =  (x,  x,  y,  ^, ),  then  one  of  the  three  production 


In  this  part  u^^  has  n-1  coordinates  while  and  u. 
n  coordinates.^  “ 


have 
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schedules  (4,0,1),  (3,1, l),  (3,0,2)  Is  optimal  for  p' .  In  order  to 
find  which  one  is  optimal  we  first  eliminate  the  infeasible  schedule(s) 
—  in  this  case  (4,0,1)  —  and  then  compute  the  cost  associated  with 
each  of  the  remaining  schedules. 

Estample  2: 

Suppose  in  example  1  that  we  let  p  be  defined  as  before  and  let 
x'  =  (0,1,0)  and  (-r*  ,x' ,yS2;' )  =  (-r>  x>  y»  £)•  (Observe  that  these 
definitions  require  us  to  interchange  the  roles  of  p  and  p'  since 
(p'-p)  is  a  unit  vector.)  Then  one  of  the  three  production  schedules 
(2,1,1),  (3,0,1)  (3,1,0)  is  optimal  for  p'.  Of  these  schedules,  only 
the  second  is  infeasible. 

Part  (b)  of  the  theorem  states  in  part  that  if  we  have  at  hand  an 
X  that  is  optimal  for  p,  if  one  of  the  upper  or  lower  limits  on 
inventory  is  reduced,  and  if  none  of  the  other  parameters  is  changed, 
then  we  can  find  an  x*  that  is  optimal  for  p'  by  comparing  the 
k(n-k)  schedules  x  +  i  <  k  <  j,  and  the  schedule  x.  This 

comparison  is  almost  as  easy  to  perform  as  that  in  part  (a),  even 
though  superficially  the  number  of  comparisons  may  seem  to  be  as  high 
as  “Ij7  “  1  (when  k  =  n/2)  rather  than  n-1  as  in  part  (a).  The 
reason  for  this  is  that  the  optimal  values  of  i  and  J  may  be  chosen 
Independently.  This  is  because  whatever  the  choice  of  i  and  J,  the 
inventory  on  hand  at  the  end  of  period  k  is  Yj^^-l  when  x  + 
is  the  production  schedule.  The  first  k  periods  and  the  last  n-k 
periods  can  therefore  be  thouf^t  of  as  two  separate  sub-problems  for 
the  purposes  of  this  computation.  Similar  remarks  apply  when  we  are 
given  an  x'  that  is  optimal  for  p' . 
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To  summarize,  theorem  3  provides  a  simple  procedure  for  finding 
an  optimal  x'  for  p'  once  we  are  given  an  optimal  x  for  p 
provided  that  +  (p-p')  is  a  unit  vector.  The  procedure  generalizes 
easily  to  situations  in  which  +  (p-p')  is  not  a  unit  vector  hut 
where  the  other  hypotheses  of  the  theorem  are  retained.  The  technique 
is  to  successively  add  and  subtract  unit  vectors  from  p  until  p' 
is  obtained.  Optimal  schedules  are  found  in  order  for  each  inter¬ 
mediate  parameter  set  using  theorem  3*  Formally,  denote  by 
p  ,p  ,...,p  the  sequence  of  intermediate  parameter  sets  produced 
in  the  process  of  successively  modifying  p  by  adding  and  subtracting 
unit  vectors.  Upon  letting  p°  =  p  and  p™  =  p',  we  see  that  each 
p^(l  <  1  <  m)  is  obtained  by  adding  to  or  subtracting  from  p^’^  a 
suitable  unit  vector.  ^ 

If  we  have  at  hand  an  x^  that  has  integral  coordinates  and  that 

is  optimal  for  p*^(=p),  we  can  determine  In  order  of  sequence 

1  2  Di  1  2  m 

x^,x‘^,,,>,x  of  optimal  production  schedules  for  p'^,p  ,  ...,p 

respectively.  Each  x  (l  <  1  <  m)  is  formed  from  x  by  applying 

the  appropriate  part  of  theorem  3. 


The  choice  of  the  p^  is  essentially  arbitrary.  Thus,  if  optimal 
production  schedules  are  needed  for  several  parameter  sets,  it  may 
be  desirable  to  define  the  sequence  Cp^}  so  as  to  Include  these 
sets.  This  may  necessitate  adding  and  subtracting  the  same  unit 
vector  In  the  course  of  the  computations.  This  creates  no  diffi¬ 
culties  but  does  Increase  the  number  of  intermediate  parameter  sets 
and  hence  the  amount  of  computation.  For  a  sequence  {p^]  in 
which  no  unit  vector  is  both  added  and  subtracted,  the  number  of 
intermediate  parameter  sets  will  be 

n 

m-1  =  2  IPi'-pJ-l  • 

1=1 
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5 


Two  questions  remain  unanswered  by  the  above  discussion.  First, 
how  can  we  find  a  p°  for  which  an  optimal  x°  is  easily  determined? 
Second,  what  can  be  done  if  we  find  an  intermediate  parameter  set  p^, 
say,  that  is  infeasible? 

We  defer  discussion  of  the  second  question  until  after  we  have 
proved  theorem  3*  One  answer  to  the  first  question  is  to  define  p° 
so  that  only  one  feasible  production  schedule  exists.  That  schedule 
is  necessarily  optimal  and  provides  a  starting  point  for  the  computa¬ 
tions.  Two  illustrations  of  this  idea  are  given  below.  Let 

_  1  _  1 

X.  1=  J  X.  and  X.  «•  V  x  . 

J=1  ^  J=1 

Start  l! 

Suppose  that  x°  >  0  is  an  arbitrary  production  schedule  with 

integral  coordinates  for  which  X°  =  R'  .  In  many  cases  it  will  be 

n  n 

natural  to  let  x°  be  a  "good  guess"  at  an  optimal  solution.  If  we 
choose  x°  =  x°  =  x°,  r°  =  r',  y°  =  max(y^,y°)  and  =  min(^|,y°) 
for  all  i,  then  x°  is  the  only  feasible  production  schedule  for 
p  =  (-r  ,  X  ,  X  ,  y  ,  y  ). 

Start  2; 

Suppose  that  x°  =  x'(x')  and  that  r°  is  such  that  r°  <  r' (r°  >  r' ) 
and  R°  =  X»(X^).  Let  x°  =  x' (x°  =  x'),  y°  =  max(y^,y°)  and 
y°  =  min(y^,y°)  for  all  i.  Again  x°  is  feasible  and  no  other 
schedules  have  this  property. 


I 


! 
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The  above  starting  points  are  intended  to  be  suggestive  and 
naturally  do  not  exhaust  all  the  possibilities.  For  example,  an 
obvious  counterpart  of  start  1  is  to  set  =  y®  =  y°. 

We  now  give  a 

Proof  of  theorem  3-part  a; 

We  consider  only  the  case  where  an  optimal  x  is  given  for  p. 

The  proof  of  the  other  case  is  similar  and  is  therefore  omitted. 

By  the  first  assertion  of  theorem  1  -  and  this  is  the  key  point 

-  there  is  an  x'  that  is  optimal  for  p'  and  that  satisfies  x*  >  z, 

where  z  =  x-(x-x* )-(x“X' )•  The  remainder  of  the  proof  consists  in 

showing  that  g(x*,r')  is  linear  in  x'  for  all  feasible  x'  for 

which  x'  >  z.  The  theorem  follows  easily  from  this  fact. 

1 

It  is  convenient  to  let  Z,  =  Y  z,  for  all  i.  We  show  that 

^  j=l  ^ 


(13) 


R'  =  Z  +  1  , 
n  n 


There  are  two  oases,  r'  >  r  and  r'  =  r.  In  the  former  event 

R'=R  +1  =  X  +1  =  Z  +1. 
n  n  n  n 

In  the  latter  case 


In  view  of  (2),  (3),  (15),  and  the  fact  that  x'  >  z,  we  have 
(ll<-)  max(x’,Zj^)  <  x^  <  min(x^,z^+l)  for  1=1,2, ...,n 
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and 


(15)  max(y^,Z^-Rp  <  X|-R^  <  min(y^,Z^-R|+l)  for  i=l,2, ...,n-l  . 

Since  p,  p' ,  and  x  have  integral  coordinates,  the  upper  and  lower 

limits  in  (ll^)  are  integers,  and  they  differ  at  most  hy  one.  A 

similar  remark  is  applicable  to  (15)*  Hence  by  P,  there  are  numbers 

a  ,a^ , . . . ,a  such  that 
0'  1'  '  n 

n 

g(x',r')  =  a  +  V  a  x' 

®  i=l  ^  ^ 

for  all  feasible  x'  >2.  Denote  by  S  the  set  of  all  indices  1 
for  which  there  is  a  feasible  x'  >  z  with  x^  >  z^.  For  these  i, 
x'  =  z  +  u^  is  also  feasible.  The  set  S  is  not  empty,  because 
otherwise  no  feasible  x*  would  exist  by  theorem  1,  contradicting  an 
hypothesis  of  theorem  3.  Now  let  1  be  an  integer  in  S  for  which 
a^  is  a  minimum.  Clearly  x*  =  z  +  u^  is  then  optimal,  which  proves 
part  (a). 

Before  proving  part  (b)  we  digress  to  establish  a  lemma  that  will 
enable  us  to  sharpen  the  results  of  theorem  1.  Consider  the  problem  of 
finding  a  vector  x*  that  minimizes  a  convex  function  f(x)  subject 
to  the  constraints 

(16)  a^x  <  bj,  1=1,2, ...,n  . 

where  the  a^  and  x  are  vectors  and  the  b^  are  scalars.  Let  S 
be  the  set  of  indices  1  for  which  a^x*  <  b^.  Also  consider  the 
same  problem  where  we  replace  the  b^  by  b^,  i.e.,  (16)  becomes 

(17)  a^^x  <  bj  ,  1=1,2,  ...,n  . 
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Lenma  1 


(a)  If  for  i^S  and  b^  >  a^x*,  for  leS,  then  x*' 

minimizes  f(x)  subject  to  (IT)* 

(b)  If  b^  <  bj  for  some  J)^S,  and  if  bj  =  b^  otherwise, 

then  there  is  an  x'  that  minimizes  f (x)  subject  to  (l?)  for  which 

a.x'  =  b*  (provided  that  (17)  is  feasible). 

U  J 

(c)  If  b^  >  bj  for  some  J)^S  and  if  b^  =  b^  otherwise, 

then  there  is  an  x’  that  minimizes  f (x)  subject  to  (l?)  for  which 

Proof; 

We  begin  with  part  (a).  Suppose  the  contrary,  that  is  x' 
satisfies  (1?)  and  f  (x* )  <f(x*).  Since  x*  satisfies  (17),  bo 
does  x"  =  0!  X*  +(l-a)x’*  for  0  <  0!  <  1.  Denote  by  B  the  set 
of  indices  i  for  which  a^^x'  >  b^.  Clearly  B  is  not  empty  for 
otherwise  x*  could  not  minimize  f(x)  subject  to  (16).  Let 


0!  =  min 
IsB 


a^(x'-x*)  ' 


Observe  that  0  <  a  <  1.  Also  for  i^B, 

a^x"  =  a  a^x’  +(l-a)ajX*  <  a  bj.+(l-Q!)b^  =  b^ 

and  for  leB 


a^x"  =  a  a^x*  +(l-0()aj^x*  =  a  a^(x'-x*)+  a^x* 
<  (bj-a^x*)  +  (a^x*)  =  b^  . 
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Therefore  x"  satisfies  (l6)  and 


f(x")  <  a  f(x')  +  (i-a)f(x*)  <  a  f(x*)+(i-a)f(x*)  =  f(x*)  j 

which  contradicts  the  optimality  of  x*. 

We  also  prove  part  (h)  by  contradiction.  Thus  suppose  that  x' 

minimizes  f(x)  subject  to  (17)>  that  a.x'  <  b'  ,  and  that 

0  J 

f(x')  <  f(x)  for  all  X  satisfying  (l?)  for  which  ajX  =  b^.  Let 

b  -b' 

x"  =  a  x'  +  (l-a)x*  where  0!  =  .  . 

Clearly  SjX"  =  b^  and  a^x"  <  b|  for  1;^  J.  Also  f(x')>f(x*). 
Therefore 

f(x")  <  a  f (x*  )+(l.Qi)f (x*)  <  f(x') 

which  Is  a  contradiction.  This  completes  the  proof.  The  proof  of 
part  (c)  is  similar  and  Is  omitted. 

As  an  example  of  the  way  in  which  part  (a)  of  the  lemma  can  be 
applied  to  our  problem,  suppose  that  the  production  schedule  x  is 
optimal  for  p  and  that,  say  y^^  ^  optimal  for 

all  for  which  y^  <  y^<  Analogous  remarks  apply  to  variations  of 
and  Xj^.  The  usefulness  of  parts  (b)  and  (c)  of  the  lemma 
will  become  apparent  in  our  proof  of  part  (b)  of  theorem  3. 

Proof  of  theorem  3-part  b: 

As  usual,  we  consider  only  the  case  where  an  optimal  x  is  given 
for  p. 

As  a  preliminary  we  show  that  there  is  an  x'  that  is  optimal 
for  p'  and  that  satisfies 
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(18) 


and 

(19)  for  1=1,2,  ...,k  . 

ae  property  (18)  follows  from  parts  (a)  and  (b)  of  lemma  1.  As  a 
consequence  If  >  yj^,  we  may  instead  let  =  y^^  without  loss 
of  optimality. 

The  proof  that  (19)  holds  for  the  revised  parameter  set  consists 
in  showing  that  for  any  x*  that  is  feasible  for  p'  (revised)  and 
for  which  (19)  does  not  hold,  there  is  an  alternative  schedule  x' 
that  has  the  properties  (i),  (il),  (ill)  (let  z  =  x  in  (li))  given 
in  the  proof  of  theorem  1.  Once  this  construction  is  Justified,  (19) 
follows  from  an  obvious  adaptation  of  the  corresponding  part  of  the 
proof  of  theorem  1. 

Suppose  x*  does  not  satisfy  (19).  Then  there  is  an  Integer 
i,  1  <  i  <  k,  for  which  x|  >  x^.  We  consider  two  cases,  yj_^  <  y^_^ 
and 

Case  1:  <  y^_^  . 

Let  J  be  the  largest  integer  smaller  than  i  for  which  x|  <  x^-. 
Such  an  integer  must  exist  for  if  not 

i-1  1-1 

y|_i  =  S  “  ^i-i  ' 

t=i  t=i 

which  is  a  contradiction. 
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Let  6  =  mln(x*-x^,  ^1-1"^  1'  ^  >  0  an<3. 

(20)  x*-x^  >  ej  i^t”^  -  ^  ^  “  jj  J+lj . .  .,1-1  and  Xj-xJ  >  e  . 


Consider  the  alternative  production  schedule  x'  defined  hy 
x'  =  X*  +  €(uj-u^).  It  is  easy  to  show  that  x'  is  feasible  for  p' 
by  using  (20)  and  the  feasibility  of  x  and  x*  respectively  for  p 
and  p' . 

Now  in  the  usual  way 


g(x»^,r)-g(x',r)  =  Oj(x|)-Cj(x*+e) 

1-1 

+  I  [h^(y|)-h^(y|+e)]  +  c^(x*)-c^(x*-6) 
t=J 

i-1 

>  €(-D"e^(Xj)  -  ^  ^ 


The  final  inequality  follows  from  the  fact  that  for  0  <  B  <  e  , 
x(8)  =  x-8(u^-u^)  is  feasible  for  p  and  g(x(8),r)  >  g(x,r). 

We  have  now  established  properties  (l)  and  (ll).  Ftroperty  (ill) 
follows  using  the  same  argument  given  in  proving  theorem  1. 


Case  2!  y|_^  >  yj_^ 


Let  J  be  the  smallest  integer  greater  than  i  and  not  exceeding 
k  for  which  x|  ^  integer  j  exists  because  If  not, 


k  k 

^  ^i-1  ^  "  ^k  ' 


t=i 


t=i 


which  contradicts  (l8). 
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Let  e  =  mln(x|-x^,  Xj-x^).  Clearly  e  >  0  and 

(21)  ^  ^  ^  t‘=i,i+l, . .  • ,  j-l>  and  Xj-x*  >  e  . 

The  alternative  schedule  x*  =  x*+'e(uj-u^)  is  feasible  for  p' 
as  can  be  shown  using  (2l).  Also 

g(x*,r)-g(x',r)  =•  Cj^(x|)-o^(x*-e) 

J-1 

+  [h^(y|)-\(y:^-e)]  +  Cj(x*)-Cj(x*+e) 

+  + 

>  e[D  c^(x^)  +  2^  D  h^(y^)-D  Cj(Xj)]  >  0  . 


The  final  inequality  follows  from  the  fact  that  for  0  <  8  <  e  , 
x(8)  =  x-8(uj-Uj)  is  feasible  for  p  and  g(x(8),r)  >  g(x,r).  This 
establlBhes  properties  (i)  and  (il)i  Rroperty  (lii)  follows  in  the 
usual  way. 

Employing  theorem  1,  lemma  1,  (18);  and  (19) >  we  see  that  there 
is  an  x'  with  the  property  that 


5  ^k  <  ^k 

and 


and 


x^-1  <  x|  <  x^  for  i=l,2,...,k  , 


x^  >  x^  for  i=k+l,...,n  . 
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is  linear  In  x'  for 


But  by  the  hypotheses  of  theorem  g(x',r) 

all  x'  satisfying  the  above  constraints,  i.e.,  there  are  numbers 

a  ,a. , . . , ,a  such  that 
0  1  n 

n 

g(x',r)  =  a  +  J  a  x'  . 

°  1=1  ^  ^ 

Let  be  the  set  of  Indices  i,  1  <  i  <  k  ,  for  which  there 
is  a  feasible  x'  with  x^  <  x^.  Let  Sg  be  the  set  of  Indices  j, 
k  <  J  <  n,  for  which  there  is  a  feasible  x'  with  ^  Since 
there  is  a  feasible  x*  either  and  Sg  are  empty,  in  which 
case  x*  =  X  is  feasible  for  p*,  or  and  Sg  are  non  empty. 

Now  if  we  fix  y^  at  any  feasible  value  between  and  yj^, 

we  may  minimize  g(x',r)  by  finding  an  Integer  i  in  for  which 
a^  is  a  maximum  and  an  integer  j  in  Sg  for  which  aj  is  a  minimum 
and  letting  x'  =  8(x',r)  =  g(x,r)+(yj^-y^)(aj.a^) . 

Thus  to  minimize  g(x',r)  we  let  y^  =  yj^,  i.e.,  x'  =  x,  if 
aj  >  a^^,  and  let  y^  =  ^k“^’  i.e.,  x'  =  x  +  if  Bj  <  a^. 

This  completes  the  proof. 

An  Algorithm  for  Automatic  Computation 

We  now  describe  a  procedure  that  starts  with  an  x  that  is  optimal 
for  p  =  (-r,x,  X,  y,  and  that  has  integral  coordinates.  The 
algorithm  then  proceeds  to  find  an  x'  that  is  optimal  for 
p'  =  (-r' ,x' ,xSySi;' )  and  that  has  Integral  coordinates,  or  discovers 
that  p'  is  not  feasible.  The  algorithm  involves  three  phases.  In 
phase  1  the  Inventory  and  production  constraints  are  relaxed,  one  at 
a  time,  until  an  optimal  x  is  found  for  the  corresponding  parameter 
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set 


p  =  (-r,  X,  X,  y,  y)  that  has  the  property  that  (x,  y)  >  Gc'>y' ) 
and  (x,  2)  <  (x',  jr')*  infeasibllltles  can  arise  In  phase  1.  The 

second  phase  Involves  moving  r  "tovard"  r' .  This  is  done  so  that  at 
all  times  each  intermediate  req.iiirements  schedule  r,  say,  is  such  that 

X  “R  <0.  This  phase  terminates  with  an  x  that  is  optimal  for 

-n  n  “ 

«  ^  M  ^  A 

^  =  (-r*,  X,  X,  y,  ^),  provided  that  p  is  feasible.  A  discussion  of 
how  infeasibllltles  are  dealt  with  in  phase  2  is  deferred  briefly.  The 
third  phase  involves  tightening  the  production  and  Inventory  constraints 
until  p'  is  reached.  If  at  any  step  in  this  final  phase,  an  infeasible 
parameter  set  is  found,  then  p'  is  not  feasible. 

We  remark  that  with  start  1,  phase  2  (and  often  phase  3)  is 
omitted.  On  the  other  hand  with  start  2,  phase  1  (and  often  phase  3) 

Is  not  needed. 

An  Important  simplification  is  possible  in  phases  1  and  3  whenever 
a  production  or  Inventory  constraint  is  not  binding.  The  idea  is  to 
take  advantage  of  part  (a)  of  lemma  1.  For  example  if  in  phase  1  we 

^  (S/  ^  ^  eg 

find  an  x  that  is  optimal  for  p  =  (-r,  x,  x,  y,  and,  say, 

y^  '  then  we  may  immediately  let  y.^  =  y^  and  still  be 

assured  that  x  is  optimal  for  the  revised  p.  Similarly,  if  we  find 
in  phase  5  an  x  that  is  optimal  for  p  =  (-r',  x,  x,  y,  £),  and,  say, 

— i  ^i'  — i  '  then  we  may  immediately  let  x^^  =  mln(x^,x|)  without 
disturbing  the  optimality  of  x.  These  shortcuts  avoid  numerous  appli¬ 
cations  of  theorem  3  and  are  well  worth  using. 

It  remains  to  develop  a  means  for  dealing  with  infeasiblllties  in 
phase  2,  In  particular  suppose  that  we  seek  an  x  that  is  optimal  for 
p  =  (-r,  X,  X,  y,  ^).  We  assume  that  p  is  feasible  and  has  integral 
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coordinates,  and  that  P  holds.  Our  approach  to  solving  this  prohlem, 

hereafter  called  problem  A,  is  to  construct  a  modified  problem,  called 

problem  B,  in  which  the  parameter  set  is  p  =  (-r,  x,  x,  y,  where 
a  a  ~ 

X,  y,  and  ^  are  arbitrary  but  not  effective  bounds.  In  problem  B 
the  cost  fimctlon  is 


where 


g(x,r)  =  X  Ci(x.)  +  ^  h.(y  ) 
i=l  ^  ^  1=1  ^  ^ 


CjUj) 


> 

^1 

Ci(Xi)  - 

)•  M(Xj^-x^)  , 

^1 

<Xi 

^i 

<Zi 

> 

hi(yi)  ^ 

'■  M(y^-y^)  , 

<^1 

and  M  is  a  large  positive  constant.  Notice  that  if  x  is  feasible 
for  p,  then  x  is  also  feasible  for  p)  furthermore,  g(x,r)  =  g(x,r). 
Let 

a  =  max  D‘c .  (x. )  , 

1  <  i  <  n  ^ 


b“  =  max  D“h  (y. ),  b^  =  max  (Zi )  > 

l<i<n  ^  ^  l<i<n 


c  =  max  _  c. (x.),  c.  =  min  _  c.(x.)  , 

X.  <  X,  <  X.  X.  <  X.  <  X. 

— i-'l—  i  — i—  1—1 


h.  =  max  h.  (y. ),  and  h  =  min  _  h  (y.)  . 

Zi<yi_<y±  Zi<yi<yi 


■  i 


1 


Let  M  Le  any  number  for  which  (h  =  h  =0) 

n  -n  ' 

+  ” 

M  >  max  {a,  b  ,  b  ,  J }  • 

i=l  ^  1-1-1 

Since  M  is  greater  than  a,  b",  b^,  the  0^  and.  h^  are  convex 
functions.  Hence,  the  modified  cost  functions  satisfy  P. 

It  follows  from  theorem  3  and  the  discussion  thereafter  that  there 
exists  an  x  that  has  integral  coordinates  and  that  is  optimal  for 
problem  B.  We  now  show  that  either  x  is  feasible  (and  hence  optimal) 
for  problem  A  or  that  problem  A  has  no  feasible  production  schedule. 

The  proof  is  by  contradiction.  Suppose  that  x  is  not  feasible 
for  problem  A,  but  that  x*  (say)  is  feasible  for  that  problem,  Then 
for  some  i,  y^  >  y^^+l,  y^  <  or  >  Xj^+1  since  the  y^, 

y^,  Xj,  and  x^  are  integers.  Thus, 

^  n  n 

g(x,r)  >  M  +  2  >  I  (°^+ii. ) 

i=l  ^  ^  i=l  ^  ^ 

>  g(x*,r)  =  g(x*,r) 

which  contradicts  the  optimality  of  x  for  problem  B.  The  proof  is 
complete. 

In  carrying  out  the  computations,  it  is  not  necessary  to  know  a 
suitable  value  of  M,  Instead,  when  a  collection  of  new  schedules  is 
examined,  as  in  applying  theorem  5,  one  proceeds  as  follows.  First 
locate  the  "least  infeasible"  schedules  in  the  class  being  examined, 
i.e.,  the  schedules  for  which  the  svira  of  the  amounts  by  which  the 
and  yj^  respectively  exceed  the  and  y^,  and  the  amount  by  which 
the  fall  below  the  is  minimal.  Then,  letting  M  --  0,  choose 
a  minimal  cost  schedule  from  among  those  located  using  g  as  the  cost 


function. 
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Ven  Maumann  Hall 
Princeton,  Now  Jariay 

Library 

Daparlmant  of  Induilrlal  Enginaring 
and  Oparatloni  Reaaarch 


Military  Mtdieal  SucMy  Agtt 
3rdAvmtMand2MiSM 
Brooklyn  32,  New  Yorit 
Attm  Mr.  Leon  Jenrtak 


Collage  of  Engineering 
New  Univerilty 
Onlveniiy  Helghti 
New  York  53,  New  York 


Librarian 

Numvlcal  Analyalt  Rtaaarch 
Unlvaralty  of  ultlornla 
Lot  Angilai  24,  California 


Logitlloi  Raiaarch  DIvltlon 
Attm  MCFR 

WHghl-Pattrion  AFB,  Ohio 


Loglitici  Raiaarch  Project 
George  Wathinglon  Unlvaralty 
707  22ndSlraal,  N.  W. 
M^lhlngton  7,  D.  C. 


Mllwadcae  Air  Procuramant  Dlit. 

AUnt  Quillty  Control  Olv. 
770  N.  Planklnton  Avonuo 
Milwaukee  3.  Wisconsin 


N.  A.  C.  A. 

1512  H.  StrHt,  N.  W. 
Waihingten  25,  D.  C. 

Attm  Olv.  of  Rn.  Information 


Newark  Air  Proeirtmnt  Dlit. 

Attm  Quality  CmM  DIvltlon 
218  Market  Slrnt 
Newark,  NawJanoy 


Now  York  Air  Proeurownt  Dlat. 

Attm  Quality  Control  DIvltlon 
niEait  16ih  Street 
Now  York  3,  NtwYork 


Office,  Ant.  Saeralary  of  Oafanta 
(RAE) 

Room  3E1068,  Tho  Pentagon 
Wnhlngion  25,  D.  C, 

Attm  Taehnioal  Library 

Operatlent  Anafyilt  Offleo 
Ntadquartori,  nolflo  Air  Porcot 
U.  S.  Air  Foroo.  APO/953 
Attm  C.  E.  ThcApion 
Senior  Analytt 
San  Franelioo,  Cailtenila 


Offloa,  Chief  of  Enginoart 
Dtpartmnt  of  the  Amy 
Withington  25,  D,  C, 

Alim  ProeiMHiint  Divlilan 
Military  Supply 

Offloa  of  Opnrailon  ANiyili 

gWXid. 

Arlln^enHall  SUtlort 
Arlington,  Virginia 
Ami  Dr.  I.  R.  Horihnar,  Jr. 


Officer  In  Cfiargt 

U.  S.  Navy  Central  Tcrpade  Offloa 
Newport,  Rhedf  Itland 
Attm  Mr.  Q.  B.  Hibloht 


Ordnanei  Carpi 
Induitrlai  Engineering  Olv. 
Dlifflond  Ordnance  Fuit  Lab, 
Wiihtngten  25,  D.  C. 

Ordnanoi  Mlialon 
While  Sands  Proving  Ground 
Las  Cruses,  Now  Maxlee 
Attm  Mr.  Paul  G.  Cox 


Philadelphia  Air  Prooifomont  Dlit. 

AUDI  Quality  Control  Olvtilon 
1411  Mlnut  Stroat 
Philadelphia  2,  Panniyivinia 

Phytloal  Rataareh  Branch 
Evans  SIcimI  Lob.,  SCEL 
Balmar,  New  Jarsay 
Attm  Mr.  Josa^  Walnslatn 


The  RAND  Corporation 
1700  Main  StrMt 
Santa  Monica,  California 
Attm  Library 


Kt  Library 

rally  of  California 
Los  Alamos  ScItnllFlo  Lab. 
P.  0.  Box  1683 
Los  Alamos ,  New  Mexico 
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Rook«l  Group 

Rodileno  Artoiutl 
HunUvItlt,  AlaUiM 
AUnt  U.  SXiBOfribva 

San  FraiMlioo  Air  ProoureimiH  Dltl« 
AUnt  Quails  Control  filvliton 
Oakland  .Arnw  Tormlnali  Bldt.  1 
Wait  Grand  I  MaritloH 
Oakland  114,  (^llfomta 


Senntori  Ordnanea  Plant 
19b  Ctdar  Avamio 
Scranton,  Ptuniylvania 
Altnt  Mr.  Cart  D.  Laraon 
Chief  inipteior 


Startdafda  Branoh 
ProeuraiMni  Divlalon 

aCS  Loglatlct*  Army 

iihtnotcn  29,  D.  C. 

Attn:  Mr.  Silaa  Wllliami,  Jr. 


Stallatloal  Laboratory 
Univarilly  of  California 
Bt^alayA,  California 


Suparlntcndsnt 

U.  S.  Maval  roitaraduato  School 
Montarty,  California 
Alh>  llvary 


TMhnIeai  Information  Offieor 
Naval  Raioareh  Laboratory 
WaihlAflton  25,  0.  C. 


Toohnioal  Operation,  Ino. 

and  C.  OrR.  0. 

Nq,  Continental  Army  Cernffland 
Fort  Mmoo,  Virginia 


U.  S«  Naval  AvIonlciPaollli/ 
Indlanapelia  16,  Indiana 
AUni  Library 


U.  $.  Naval  Cnfllmerlnq  Expeilment 
Station 

Anapolli,  Mainland 
AUni  Mr.  p,  R.  t'llPrIare 


U.  9.  Naval  Inspeotoi  of  Ordnanoe 
400  S.  BelflerSveet 
Mlihawaka, Indiana 


U,  S.  Naval  Inipeotor  of  Ordnanoe 
Eattman  Kodak  Company 
90  W.  Main  Straet 
Rocheiltr  14,  New  York 


Weitern  Conk'aot  Manaeement  Reqton 
Attni  Directorate,  Qtmllty  Control 
Mira  Loma  AP  Station,  California 


Dr.  Adam  Abritail 

Dept,  of  Eoonomloa  and  Enolnooring 
Stwenc  (nitltuta  of  Technology 
Hoboken,  New  Jersey 

Dr.  Max  Aitraehan 
Dapartmant  of  Loglattc 
The  RAND  Corperttlen 
1700  Main  Street 
Santa  Monica,  California 


Profaiaor  T.  W.  Andarson 
Dopartffltnt  of  Maihefflaileal  Sutlitlci 
Columbia  Untverilty 
Ntw  York  27,  New  York 


Profeitor  Fred  C.  Andrews 
Mathematics  Department 
University  of  Oregon 
Eugene,  Oregon 


Professor  Robert  Beehhofer 
Sibley  School  of  Mech.  EnglncHlng 
Cornell  University 
(thace.  New  York 


Profeiice  R.  E.  Beckwith 

Graduate  Sohaol  of  Builneit  Administration 

UnivertUy  of  Southern  California 

Los  Angeles  7,  Callfemia  1 


Dr.  Franc li  Dretch 
Stanford  RHCiKh  Instltula 
939  Ravonswood  Avenue 
Menlo  Part,  California 


Professor  J.  N.  Beirattent 
Weiteni  Reserve  L’nlverilly 
ClevMand,  Ohio 


Mr,  Carlton  M.  Beyer 
Offlee  of  Gutdod  Miitlles 
Offke  of  Atil.  Secrotaiy  of  Defonio 
(RiEl 

WashlngioA,  D.  C. 


Professor  Arohtson  J.  Duncan 
Diparlffltm  of  Induslrlaf  Engineering 
Johns  Hopkicts  UnIvortUy 
Baltimore  18,  Maryland 

Profeitor  Meyer  Dwaae 
Department  of  Mathemalloi 
Northwestern  Univsrtity 
Evanston,  Itllnoli 


Prefetior  2.  W«  Blrnbaum 
Laboratory  of  Sutlatloal  Roiiaroh 
Department  af  Mathenatloi 
Uni^lty  of  Washington 
Seattle  9,  Waihinglen 


Prefesser  RmioII  Bradt 
Oepartment  of  Malhamatlos 
Univeraity  of  Kansai 
LawrmMO,  Kanias 


PMif«i««  irviM  W.  Burr 
Oepirtmmtt  of  Malhamatloi 
Purdue  University 
Lafayette,  Indiana 


Mr.  Q.  Burrows 
Knolli  Atomto  Power  Lab. 
Sohentotady,  New  York 


Or.  A I  Charaes 
The  Teohndoeieai  Inilltuie 
Nerthwetlern  university 
Evaniten,  iitMoia 


Mr.  W.  H.  Ciatwerthy 
Bettis  Pien,  WHtlneMuse 
Eleotfto  Cerpeeatfon 
Dm  iitt 

Pittsburgh  90,  Pennsylvania 


Professor  0.  A.  S.  Fraisr 
Oepartmsflt  of  Mathonutles 
Univariiw  of  Toronto 
Toronto  5,  Canada 


Mr.  Wllllim  E.  Ollbart,  Chlaf 
Mathematical  Statlitlos  Branoh 
Atonlo  Energy  Commlialon 
IMihIngton  25,  0.  C. 

Mr.  Lean  Gilford 
Operations  Retsaroh  Ino. 

8o09  Cameron  Street 
Silver  Springs,  Maiyland 


Mr.  uernard  P.  Goldsmith 
Asieolale  Prefeiior 
Norlhiailern  Univsrilty 
960  Huntington  Avtnui 
Boston  15,  MaiiaohuitlU 


Prefosscr  Leo  A.  Goodman 
Statistical  Rf  joaroh  Centar 
Unlvoriltv  of  Cnioago 
Chleago97,  llllneli 


J.  Onranwoed 
490  Great  Palls  Stroot 
Falls  Church,  Virginia 


Professor  Paul  Clifford 

New  Jirny  Slate  TeMhere  College 

Menteiiire,  New  Jersey 


Professor  W.  Q.  CMinn 
Oepartment  of  Stetietloi 
Harverd  University 
Cambridge,  Miisaohutilts 


Profeuor  C.  C.  Cookertwn 
Institute  oTSUItslIos 
Suu  Celleat  Section 
Nath  CeroilM  Slate  Collega 
Raleigh,  North  Car^tu 


Profess*  Edward  P.  Coieman 
Engineering  Dipartment 
Univeratty  of  ^llfornla 
Lot  Angeles  24,  California 


Miss  QartrudeM.  Con 
laitltuu  of  Suilitles 
North  Cardlna  State  College 
fUlelgh,  North  Cvellaa 


Dr.  Joeeph  Daly 
U.  $.  Censut  Bureau 
Washington  25,  D.  C. 


Professor  Cyrus  Dermsn 
Department  of  Induitrlil  Engineering 
Cdumble  Untvarstty 
New  York  27,  New  York 


Mr.  H.  F.  Dodge 
R'tgers  Universitv 
Ntw  Brimswlek,  New  Jirsey 


Professor  Donald  Caver 
Wisllnehouit  Reseereh  Labs. 
Euteh  Road 
Churchill  Bore 

Pittsburgh  99,  Pinniyivenie 


Profess*  Frank  M.  Gryna,  Jr. 
UnlveriKv  CoHvoe 
Rutgers  University 
New  Brunswick,  New  Jersey 


Dr.  Donald  Cuthrli 
Stanford  Research  Inatllute 
999  Raveniwood  Avamw 
Menlo  Perk,  Celll*nla 


Dr.  Thtod*e  E.  Harris 
The  RAND  C*ponitlen 
1700  Main  SUeel 
Santa  Monica,  California 


Dr.  Loon  H.  Harbaek 
Dapt.  of  Industrial  Engr,  &  Oparatloni 
Riiearch 

Collage  of  Engirwwing 
Now  Y*k  Unlvorslly 
New  Y*k  59,  New  York 


Profess*  W.  HIrseh 

Institute  of  Mathematical  Sciences 

New  Y*k  Unlverilty 

New  Y«k  3,  NewYwk 


Mr.  Eugene  Hixson 
Code  600.1 
GSFC, NASA 
Grecnbelt,  Maryland 
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Dr.  Paul  0.  HmI 
DfPirtawnl  af  Mathcmatlca 
Uniwtlty  of  CalMa. 
Ln  Aaimh  CalllomUi 


Profmor  Harold  HoUHInt 
Aiioolato  Dlroolor 
Initiwt  of  Sutinioi 
Unlvanity  or  North  (Sirolliu 
Chapol  Hill,  North  CarollAi 


Profniff  L.  Hinrifli 
Soh^  M  Buolnaii  Admialitrailon 
Unlvanity  of  MImoioU 
MlAAoaoolii,  Mlmaiatt 


Dr.  JiAM  Ik.  JaakiM 
MaoiMMiit  Sclaneai  Raioareh  PrejNl 
6S  AMhliVaklM  klidlAfl 
Unlvanity  a?  ffiamU 
Lot  Aiiflrtat  24,  £Ul?ornla 


Or.  W.  C.  Jaoob 

Em  OoMrtmrtt 
ilty  of  llllnolt 

I  nilfwia 


ProfonorW.  D.donto 
OoMMionl  of  Itatlitloi 
Ml^loan  SkaU  Unlvtfilty 
Baak  Lamina,  Mlehlfu 


Ccniriiif^lilon 
Qofloral  lorvlOM  Adailnlalrotlon 
Roan  b316i  Raalan  S  Building 
WtahUigton  »,I).  C. 

ProfMior  Oaear  Kofaptbama 

•liailtUot  LikopKoy 

Iowa  suit  Collnoo 
Anaa,  Iowa 

Snaor  SototMt  Kullbaok 
irtnonk  of  Stallit'o^ 

■go  Mihlnglon  U>)vonlty 
Wathlntfon  7,  u.  C. 


Or,  Carl  F.  Koiuek 
IBM  Corporation 
2401  So.  Main  Slnot 
Howlon  2,  Toxii 


Hward  Ultin 
170S  W,  Artoala 
Torranot,  California 


Or.  E.  L.  LfClarg,  Chiof 
Blomrtrleil  Sorvloai 
Agrifltilltral  Rttovoh  Sarvlet,  USOA 
Boliivllia,  Maryland 


Profotior  Sabatllan  6,  lltUuar 
411  Enginoarirtg  Building 
Columbia  Unlvanity 
Now  Voii  27 ,  Now  York 


Dr.  Eugona  Lukaca 
Dapartmont  of  Mathamitloi 
CalNia  Unlvonlty 
Waahtngton  17,  D.  C. 


Or.  Robort  Lundaoird,  Hoad 
Lodtllaa  and  MatAomaueal 
BrtMh 


imaueal  Sutlitici 


Offleo  of  Naval  RaMoroh 
Mahingion  29, 0.  C. 


Dr.  Clifford  J.  Malonm 
DIvltlon  of  Blologlea  siandarda 
Natlomt  inatituti  of  Hoalth 
blhaada  14,  Maryland 


Profottor  Frinli  MaiMV 
School  of  Publlo  Haaltn 
Univaritly  of  Caltrornta 
Loa  Angaiaa  24,  California 


Profotiar  Q.  W.  UaBlitth 
DaiartwonI  of  Induotrlal  EnglMorlns 
Unlvortlky  if  Mlamola 
Hlmaa^li  14,  MbMHOla 


Or.  Paul  Moyar 
paportmaitt  ar  hiaihMiatloo 
mihiAfton  8uto  Unlvaniiy 
Pullman,  WuhJjigkon 

Dian  Paul  E,  Mafui 
School  of  EnglMarlnQ 
Unlyarally  ofluhalo 
Buffalo,  Now  York 

Mr.  R.0.  Mwghy 
Bali  tilaaheno  Labaraurloi, 
4bS  Wairstraal 
Nm  Vtrt  14,  IMw  Vark 


0  B  Novndiam 

UI.I/  MmIiI.1  EiiluHrlii.  BIvIiIm 
CaoMrallof 

Hdgn. ,  San  Bamaidino  Air  Mataalol  An 
Nanart  Air  Faroo  Bait,  Callicmia 


ProfHiv  J.  Noyman 
Barfcaiiy4,  CaMfamla 
Mr.  Prad  Okano 

Natianal  AaraMutiai  k  Bpooa  AdninliiraUan 
RallSiiTll^  Sytloma AnaiyiltOITIaa 

SA-lSl 

I  OH*  Straotj  N.  W. 

Ingion  ft,  D.  6. 

Proftiiar  B.  R.  Ott 
PMArt— lU  of  lulhNWtlfit 
Rulgm  Unlvorallv 
Now  Brunpwlok,  Now  Jaaopy 

Mr.  Cyril  Poekkam 
PN^lOlabt 
Untvan)^  of  Dayton 
Dayton,  Ohio 


Or.  Richard  Poit 
Dtpartmtnl  of  Mathamillot 

?anJoit  Slat#  Collagt 
an.)oat,  Calllemta 

Prefaiiar  P,  H.  Randolph 
Purdua  Unlvinlty 

Dtparlfflont  of  Inouilrlal  Englntarlng 
LafaytUa,  Indiana 


Prolfiier  Otorgo  J,  Rtanlkoff 
DMartmant  of  InduaUlal  Englnoorlng 
lliimli  Inalltuta  of  Tiohnoiogy 
Chisopo  16,  ilMnola 


Dr.  Paul  R.  RWor 
Chiaf  Statiatieian 

Aaroniullaa)  Roitarch  Lab.,  WADC 
WHgM-Paltaraon  AFB,  Ohio 


ProfMior  Harbart  Robblm 
Malh^leal  Stattillai  Dtpartmanl 
Ceiiffflbli  UffIvcriHy 
NowYwk  27,  Now  York 


Or.  Harry  6.  Romlg 
391  Aira  RmI  Drlrt 
Paolflo  Palloadai,  Callfomta 


Dr.  Harry  Roionblatt 
SUtlillaal  Rtiaaroh  DIvltlon 
Bureau  of  Camui 
Waahln^n  29,  D.  C. 


Manor  HtTMA  Rubin 
Oopwtmont  of  StatlMloi 

ProfMior  NinM  Rudy 

BUtlillei  OMorbwHil 
SaaruiMto  Siau  Colloga 
Saoramanto,  California 


Proftiiar  Murray  Roifnblatt 
Dopartiutnl  of  Mathtmtlei 
Brown  Univoiilty 
Providonea,  Rhode  Itlend 


Mill  Marlon  M.  Snndomlri 
2ftl  Cedar 
Baneley  9,  CallfamU 

Mmiot  I,  RIehard  ^vaie 
Sohool  If  BuilnMi  Aemlnrun 
Unlvanity  of  MliwMon 
MInnaapMli,  MlnnMOlA 

Profaiiv  L.  J.  Saviga 
Mathaiwiloi  Deportmonl 
Unlvanity  of  Mkhlgin 
Ann  Arbor,  Miehtgan 

ProfMior  Hanry  Bohtff^ 
OaMMiM  of  Sttllfiioi 
Unlmlky  of  CMifomla 
BirkaiiyA,  Garifomla 


^ofMior  Robort  Bahlalfer 
Giiduato  Sohool  or  BmIahi 
AdmlnliiratloA 
Harvard  Unlvanity 
BmIoa  63,  MaiiaahuMUa 


BrlgadlvGoMrai  W.F.B,  Sahradar 
e^U.  Cal.S.t.C.eiMawli 
Mloa  if  U«  Military  AHaoha 
Auftnilan  Bmbaiiy 
200IComMlioiit  Am.,  N.  W. 
Bm4837  ^ 

WiaMffglonBfD.  C, 


Profnaor  Baynieur  ShwiM 

StmoiH  of  Mathomatloi 
I  Stole  Univmity 
It  2,  MIohlgin 


Mr.  Wiltar  Shmuhirt 

Bell  Ttitphono  LabontoriM,  Ino. 

Mimy  HHI,  Now  Jori^ 

Dr.  Roiadllh  SItraavti 
Toaelwn  Callagi 
CMumbla  Unlvanity 
Now  York  27,  Now  York 


Dr.  Milton  Sobel 
Statlitloi  DfMrtmant 
Unlvaniiy  of  MIimmoU 
Mlnniapelli,  Mlnnaiau 


ProfMior  Frank  Spltear 
Da  partimnt  or  Mauwmitlei 
Comall  Unlvcnlty 
Ithaea,  Naw  York 


Mr.  Arthur  sum 
Comil  Aaroniulloal  Lab.,  Ino. 
P.  0.  Ban  235 
BufU}o»,  NowYvk 


ProfMior  W.  Allan  Willli 
Offloa  of  tiw  Pnildtnt 
Unlvanity  of  Roohoitar 
Roahaitar  20,  Now  Y«k 


Dr.  Harry  Walng«Un 
United  StaUi  Armi  Control  end 
OliariMment  Agency 
Waihington  ft,  D.  0. 


Or.  Irving  Wclai 
The  MItn  Cerpcntlon 
Bedford,  MaiiaehuseUt 
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1 


ProTeitor  J.  S.  Rustaal 

Cell«9«  of  M^lelnt 

Univarslty  of  Clnehnall 

Clnolimatt,  Ohio 

1 

Or.  M.  0.  Ktfldall 

C.EJ.R.,  UKoLtd. 
furrliA  BviMlai 

Groat  kVHt  ItMi 

BrffMfoid/MMda. 

1 

Mhm,  VMoata 

Capuin  Burton  L.  Willar 

AF^  Offlea 

Martin  Aircraft  Corporation 

Danvor,  Colorado 

1 

DPMtar,  iMtRai  da  SimIsUrm 

9,  QmI  UM  Bamaid 

Paris  V,  t^ranoa 

1 

I 

Profatior  Mason  Wasoott 

Editor,  Industrliit  Quality  Control 
Riitaars  Univortrtv 

Now  Brunswlok,  Now  Jmay 

1 

Dr.  ftudoft*  Hussar 

lastttut  f»  Anaowandta  Mithamilli 

IMIvmlUt  Bam 

HoehsaluiltiraBsa 

Bant,  Swhiariand 

1 

Mmir  0.  S, 
y>Mnlt»  VTmm. 

][gj«{|«  S,  mM.  ^ 

fvofHsor  S.  S.  Wilks 

Dapartmont  of  Mathomatles 

Prlneaton  Unlvorsily 
mnetlon,  NowJorsoy 

L 

Profassoi  Amllo  Heyland 

Oaaartiaant  of  Mathaaaltos 

Unl^lty  af  Oslo 

Blindsm,  Norway 

1 

Profoisor  J.  Wotfawiti 

Dapartmnl  of  Mathoaulica 

Comall  Univortity 

Ithaca,  Maw  York 

1 

Profoisor  A.  Willhar 

Ttohatseha  Nochtohala 

Oanaitadt,  GariMny 

1 

Dr.  N.  C.  H«tmktr 
PhIMpi  Rittirth  Labontcrlet 
Eid«nhov«n,  Ntth«rlandf 


Mr.  I.  0.  Hill 
SUtlktloal  Advisory  Unit 
Mlnitliy  of  Aviation 
Shtll*Mox  Houia 


Strand,  London,  W.  C.  2 


P,  V.  Slvaranlirlslian 
SQC  Unit,  Indian  Stat'.illoal  liwUtuU 
Qiiatni  Read,  Oovammont  OfflON 
klldlnt 
Bombay,  iMla 


SivinBaabnuihana 

SQU  tMl  Bannlora 

c/e  SKSJ  7oehftole0lcai  Imt. 


6an«aIori  1,  India 


Caniraet  Moflr«32B(93) 
StptMiMr  1962 


Vi 


PfOfiltorG,  E.  P,  6m 
DturUnefit  of  Stotittiei 
Univirilly  oTWItoMsIn 
M*dlion«  Witeoniln 


Profoitor  lUlDhA.  Bradloy 
Dopartmont  of  Statliflot 
Florida  SUU  UoiyMtHy 
TalfihMSMi  Florida 


PrefoitorW.  T.  Fodorar 
Comoll  Univorelty 
OoparUntnt  of  Plant  Groodtna 
Blomatrlot  Uolt 
lthaoa«  NoyvYork 


Prefoitor  H,  P.  Good# 
Dopartnofit  of  Induttrlal  and 
EnolnaorlnB  AdmlnlitraUon 
Camoll  Univartity 
Ithaca,  Now  York 


Profoiacr  Paul  Malar 
Dapartmant  of  Statlilloa 
Unlvaraliy  of  Chleago 
Chleaao,  llttiiota 


Profaiiop  W.  L«  Smith 
Slatlitlai  Dapartmant 
Unlvaralty  of  North  Card  ina 
Chapal  Hill,  North  Ctfdlna 


Profaiiw  M.  B.  Wllk 
Statlatloa  Cantor 
Rutgan  <•  Tha  Siata  UnivanIty 
Naw  Gruntwloh,  Haw  Jariay 


Mr.  William  Mllna 
Oaparimant  of 

Inipaailon  i _ 

7S  St.  Patrlflk  Siraat 
Ottawa,  Ontario 
Canada 


}f  National 
n  Sarvloaa 


Oafania 


Mr.  C.  F.  Hobbi 
Syatamt  Analyila  Saotlor 
Arm  Davalepmant  filabllihmant 
2S9  Argyla  Avanuo 
Ottawa,  Ontario 
Canada 


Profaatto*  W.  r.  FrallMrair 
Malhamatlot  Oaprirtmart 
Brown  Unlv«*iHy 
Provldanoa,  Rhode  tiland 


Prefaiaer  Robert  Thrall 
wnt  of  Mathamatlot 


Oapaiifflint  of  [ _ 

Unl^flty  of  Michigan 


Ann  Arbor,  Mlohlgin 


Or.  Harman  Qolditina 
IBM  Raaaaro 
Boa  218 

Yorirtown  Halgfiti,  Naw  York 


Dr.  F.  J.  Murray 
U.  $.  Amra  Raioareh 
Boa  CM,  Oidta  Station 
Durham,  North  Carol Irta 


Profaiaor  Jack  Klafar 
Dapartmant  of  Mathamitloi 
Cornatl  Unlvaralty 
Ithaoa,  Naw  York 


Or.  John  Tliehandorf 
Allentown  Laboratory 
Ball  Talaphena  Labontoriai,  Ino. 
Allentown,  Pamiylvanli 


Mr.  Martin  B.  Brouda 
Quality  Control  DIvlilon 
Ganaral  Sarvleas  Admlnlitntlon,  FSS 
Room  6919,  Reoton  9  Bldg. 
Waihingten  25,  D.  C. 


Mlai  Tharaia  Friaka 
Raiaarth  lirttltuta 
Uahranlty  of  Dayton 
Dayton  9,  Ohio 


Profaiacr  J.  V.  Hanlaon 
Qiallty  Control  Ceuraa  DIraatv 
SahM  af  L^allaa 
AroaA,  P.  0.  Box  • 
Wrlght-Patlaraon  AFB,  Ohio 


Mr.  Pater  K.  Kiiffnar 
Tha  Halhenftara  Company 
9th  A  Kettnar  Avammt 
Chicago  24,  lllinoii 


Dean  E.  A.  Trabant 
Sohod  of  Eim  ' 


Sohod  of  Eiwingaring 
Urtlvanlty  ^iMrilo 
Buffalo,  N^  York 


^pt.  W.  H.  Koon.  USM 
Oapartmaiit  of  tho  fuvy 
Ornaa  H  Naval  Riaaareh/Coda  400 
Waihingtan  25,  D.  C. 


ipMirMndlng 


OfflM  of  Nival  ftiiairoh  Brinoh  Ofnei 
lOlO  Sail  Oraon  Siroat 
Paiadana  1,  California 
Atini  Dr.  A.  R.  Liirfar 


Maihamallci  Branch 
Atmy  RHOWOh  Offlaa,  OCRD 
Papartmant  af  tha  Amy 
mihtrtglM  29,  D»  C. 


Chlal 

U.8.  Army  RAD  Llalion  Oroup 
APO  797 

Naw  York,  Naw  York 
AUni  Solantifle  Advlaor 


Librarian 

Array  Loglitloa  Managamant  Cantir 
Fart  LHi  Virginia 


Commanding  Ganaral 
U.S.  Army  Combat  Davatopmant 
Exparimtril  Cantor 
FartOrd,  Callforali 
Attni  Stanford  Raiaaroh  IniUtuta 


TMhnIeil  Library 


U.  S.  Mllltara  Aoadany 
%«t  Point,  Raw  York 


Combat  Oparat..  .  . . 

Hg.  Contlnantal  Army  Cemiund 
'**■'  ** -  *'*"1011 


Fart  Motwot,  VIrgInl 


Santer  Solanllii 
Humiit  Raiiarah  Unit  No.  1 
Contliwntal  Army  Command 
Fart  Knox,  Kentucky 


Director 

Rastireh  Anilyili  Corp. 
6999  Arlington  Road 
Bathaada,  Im. 


Human  Raiouroai  Ractireh  Offlea 
Gaerga  Wiihington  Univarilty 
Box  9996 
Wathinglen  7,  D.  C. 


Comnardlng  Qenaral 

U.  S.  Army  Ordnanea  Ammunition  ComMnd 
Jdlot,  llllnolt 

Attnt  Technical  Library 


Cotmiandlng  Oanaral 
U.  $.  Anv  Ordnanco  Special 
Waapana  Ammunltlm  CBmmwd 
Dover,  Naw  dtnay 


Commindlitg  Oanaral 


U.  S.  A^  Ordnanea  1 
Autoiietlva  Coamiand 
1901  Botfd  Siraat 


Oalrolt9,  MMIgiit 


Comniindlng  Ganaral, 

U.  S*  Army  Ordnanea  Woepang 


C.  Q,  Rook  liland  Araanal 
Rookliiand,  llllnolt 
AUnt  Ttohnioal  Library 


Diraetor 

nBilliitIo  Mliilla  Aganoy 
.  Arriv  Ordnanea  Mitillo  ' 
Raditona  Arianil,  Alabama 


OIraotar 

Am  Rookat  A  Guided  MItilla  Aganoy 
U.  9.  Araw  Ordnanaa  Mliilla  Commuid 
Raditona  Mtanal.  Alabama 
Attrii  Chlal,  RAD  DIvlilon  1 


Diraater 

I  Rocket  A  Guided  Mliilla  Agiwy 


Arm  Rocket  A  Guided  Mliillr 
U.  S.  Arnw  Ordnanoi  Mlaalta 
Raditona  Araanali  AtabaiM 


oiMna  rwitmii  mtaBamt 

AUnt  Enginaarfng  RadulraminU 


. .Raoww 

CecfdiMtian  OFF 


Diraater 

FiltMA  Raiaarah  A  BnufMorlni  Libi. 
Ploatinrw  Ananal 

Dover,  New  Jeniy  1 


DlTMter 

Mijer  Item  Supply  MiMgtmnt  Afiney 
Lattarkonny  OranaAcc  papet 
Chatnbariburg,  Panniylvinla  I 


Director 

U.  S.  Army  Orriimneo  Supply 
Analyiii  Aotnoy 
Maluehtn,  new  Janay 


Chief,  StoroM  A  Mtlnlananea  DIv. 
Dfraoierata  for  Quality  Centrd 
MGSA 

Richmond,  Virginia 


Commanding  Offloar 
Htadwirtrti  Ordnanea  Million 
WillaSindt  Mltcllt  Range 
Lai  Cruoai,  Naw  Mexloo 
AUnt  ORDES-OM-TIO-TL 


Camnandlng  Officer 
Whlarv  let  Arttnal 
^tirvilat,  Naw  Vark 


Conrnndlng  Offloer 
Radford  Ananal 
Radford,  Virginia 

AUnt  Chlifj  IntpaetlottOlv,,  ORDCB 


DIrHtor 

U.  S.  Army  Corpa  Of  Engincan 
Engtnacr  Rtiairch  and 
Devalopmcnt  Laboratories 
Fort  Balvoir,  Virginia 


Offlea  of  tha  Quartarmutiar  Oanaral 
Operational  Mithamiilci  Branch 
mihlngton  25,  D.  C. 


Contract  Nonr-BBStSB) 
Saptamlrar  1962 


I 


Commanding  Genaral 

U.  S.  Army  Transporutlon  Material  Command 
209i  Main  Otflee 

St.  UuU^  Mitaourl  1 


CtfiRtaiidfitg  OfTicor 
dU  Rotenreh  &  Engineering 
Fl^i  evaluation  Agency 
n<  L«««  Virginia 


EiiBtrenls  Syitsmt  Diviilon 
Atuii  Quality  Cofi^ 

L.  Q.  Hanacoffl  Field 
Bedford,  MasMohuieUt 


Technical  Library 

U.  S.  Army  Electronic  Proving  Ground 
Fort  Huachuca,  ^rlxena 

Teehnloal  Library 
U.  S.  Army  Sigul  Refoareh  and 
OtvelmMot  Laboratorlei 
Port  Monmouth,  New  Jersey 


Oiltf,  teohntctl  Evaluation  Braneh 
Aiv  Reaotiea  Office,  ORB 
(1,  S.  Alonta  Eftergy  Coemlitien 
tfahlngton  2S,  D.  C. 


iliihomitlee  Rtteanh  Canter 
tlnVentiy  cf  WIscenaIn 
Uafiaen,  Wlaoentln 


Rome  Air  Pone  Depot 
Management  RequIremenU  Dhr. 
Management  Analyili  Section 
Qrlfflii  Air  Force  But,  Mewyorfc 


Rome  Air  Force  Depot 
(klfflit  Air  Foreo  fine,  I 
Attni  W 


Conunandlng  Officer 

U.  $,  Army  Signal  Clectrcnlca  Research 
Unit 

P.  0.  Box  20S 

Mountain  View,  Califerau  1 


Commanding  Qefwral 
U,  S.  Army  Chemical  Center  Corps 
Engineering  Command 
Army  Chtmloal  Center,  Maryland 
Attni  Mr.  T*  M,  Vlnlng 


Director 

U.  S.  Army  Chemioai  Resoanh  and 
DovelqprMnt  Laboratories 
Arm  ChemlenI  Contor,  Maryland 
Aunt  TeohnieaJ  Library 


CooMiandlng  Olfloor 

U.  S.  Army  Chemlul  Corps  Proving  Ground 
Dufway,  Urah 

Coom^lni  ONloar 

New  Torh  Chemioai  PreeieoeMnI  DIat. 

UO  Varlok  Street 
NewVorkld.  Hew  York 
Atini  Quality  Surely  Division 

MatherMtlus  Ulvlilen 

Air  Faroe  Office  or  Seientifle  Resuaroh 

Withington  29,0.  C. 


OirMlev 

Censdian  Army  Ooirallonal  Rotooroh  CsUb. 
Oepartment  of  National  Defense 
Otuwi,  Ontario.  Ciraidi  1 

Attni  6.  A.  Rtclardion 
C.  A.  Do  Cande  e  CA  ORC/8/2900-1979/S 


Statiitlos  Department 
George  ttashington  Uni 


Personnel  Research  Branch 
Personnel  Reiesrch  li  FVooedtfi 
The  Adiurant  Generals  Offloe 
Wishbtgton  29,  0.  C.  • 


OIreolor 

Jet  Prg^slon  Laboratory 
4600  Oak  wove  Driva 
Pasadena  9,  California 


National  Inslltuto  d  Hosllh 
U.  S.  Publle  Health  Servloe 
Bolhnda,  Maryland 


Utrarlan 

ill  S.  A.  Artlllary  A  Mlitlla  Sehool 
Pel  Sill,  Okialiena 


Cttnrandlnt  OeMfil 

Cramand  and  Oeneral  Staff  Collags 

PI.  Lmavonwoftli,  Kaniai 


CMmanding  Offleir 
Pkitlimy  Araanal  /O 


CramandlM  Offleor 

(manes  MaUrlali  RHoaroh  Offke 


W^erra  Arsenal 
ybtirtewn  72,  Mai 


CNMiandlni  Offissr 
Lititt  Inipeellen  Dlvlslert 


CMffMtndIns  Off)e< 
^0  Teoinlo^  ! 


(m  Llalien  Qreup,  Projeet  h 
wITow  Run  Researeh  Center 
YpKarPJ,  Mlshlgu 

fiirroter 

Llmoln  Laboratory 
Uungtef.,  Miiiaehusetti 

/tNlied  Mathematloi  Labs, 
hillofral  Bureau  of  Standards 
t^hln^on  29,  D.  C. 


TKhnIoal  Informitlon  Servlet 
Mini  Refererwo  Branch 
P.O.  BoKb2 
Ott  Hidgi,  Ttnneaite 


Oitf,  Bureau  ol  Naval  Weapons 
diMlv  f  strol  Division  /tQC/ 
Ctfartfno(.(  of  the  Navy 
V^UrgtoR  29,  0.  C. 


Cstgraphy  Depertnienl 
%  George  WaihliMton  Inlvarslly 
WHhlngtwb,  D.  T, 

AUrtt  Professor  R.  D.  Campbell 

SapsTlntendonl 

V.S.  Naval  PestgraduaU  School 
IlffUerey,  California 
Attns  Library 


Romo  Alt  Foreo  Depot 
Griffisi  Air  Force  Base,  f 
Attni  ROPQ 


Dayton  Air  Foreo  Depot 
Dayton  20,  Ohio 
Attni  MDQ 

Dayton  Air  Forte  Base 
Dayton  20,  Ohio 
AUni  m6mQ 


Daytan  Ajr  Force  Depot 
Dayton  20,  Ohio 
AUni  MDSQ 


Aeronaulleal  Systemi  Division 
Alln<  MKy  Control 
WtrlghtoPalterson  AFB,  Ohio 


2709  Air  Faroe  Vehicle  CoMrel  Group 
AtUii  Quillty  ConM  Div. 

Meii#li,  Tennesiee 


Dgden  Atr  Preewemant  Diilrlot 
Altni  Duality  Control  Olv. 

Hill  AN  Feree  Base,  Utah 

Riohmond  Air  l^urtmenl  Dlitrlot 
Attn;  Qi^Ry  DIv. 

P«  0.  BonBSoS 
Rletu^  29,  Virginia 


U.  S,  N.  0.  T.  S. 
Pasadeni  Anntii 
1202  E.  Fo#lll  Mvd. 
Patadena  B,  Callfomla 


Chief,  Bureau  ol  Yards  A  Dooks 
Depertment  of  ihs  Navy 
Mhshinglm  29,  6,  C. 

Alint  C^e  A69D 

Commanding  Offlear 
U.  9.  Naval  Avlenles  Facility 
21sl  St.  A  AHinflan  Avenue 
Indlanapolli,  Indiana 


Bureau  of  Naval  Whapent 
Branch  Rtprasantatlve 
Eaatmin  t^k  Co, 

121  LIneoln  Avenue 
Roohesler  11,  New  York 


Headouartirs 
AlrPorce  Syttsms  Commuid  /SC2KQ/ 
VWIght'PaUeroon  Air  Force  Base 
Dayton,  Ohio 


frtwam  Dirteter 
UdnarniUMi  Soltneti  Program 
Mltonat  Sotanoo  Foundation 
yhihmgten29,  D.  C. 


Dlrsctor 

Wsltor  Read  Army  Initllutt  of  Researeh 
Wkllir  Reed  Army  Medioal  Center 
IMthtngten  29,  6.  C. 


SotenUlle  OIrtoter 

9M  Resewoh  A  Engineering  Convnand 

.  8. 

Nallok,  Massachusetts 


Commandant 

AM  Food  &  Container  Inst  for  Armed  Forets 
U.  S.  Army 

1819  Pershing  Road 


I*.  r 

Chicago  9,  Illinois 


CmoHMd'ng  Offioor 

U.  S .  Army  Transportation  Rosaarch 

Cmnond 

Pot  Eustls,  Virginia 
Mtrr:  Rosearoh  Relarenee  Cantor 


HMdqMutors,  United  Stales  Air  Force 
AUnt  AFCOA 
Mihlngton  25,  D.  C. 


(Xindo  Air  Prooietment  District 
Attnt  Quality  Conl/ol  DIv. 
Cdindo  Air  Force  Base,  Florida 


Additional  copies  f«  project  leaden 
and  aisistantf ,  offlne  fife,  and 
raserva  for  futirt  roquramenU 


Contract  Nonr-225(53) 
SepteiMwr  19b2 
(495) 


